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ABSTRACT 

We present simultaneous Plancl<, Swift, Fermi, and ground-based data for 105 blazars belonging to three samples with flux limits in the soft X-ray, 
hard X-ray, and y-ray bands, with additional 5 GHz flux-density limits to ensure a good probability of a Plancic detection. We compare our results 
to those of a companion paper presenting simultaneous Planck and multi-frequency observations of 104 radio-loud northern active galactic nuclei 
selected at radio frequencies. While we confirm several previous results, our unique data set allows us to demonstrate that the selection method 
strongly influences the results, producing biases that cannot be ignored. Almost all the BL Lac objects have been detected by the Fermi Large Area 
Telescope (LAT), whereas 30% to 40% of the flat-spectrum radio quasars (FSRQs) in the radio, soft X-ray, and hard X-ray selected samples are 
still below the y-ray detection limit even after integrating 27 months of Fermi-hAT data. The radio to sub-millimetre spectral slope of blazars is 

quite flat, with (a) ~ up to about 70 GHz, above which it steepens to (a) 0.65. The BL Lacs have significantly flatter spectra than FSRQs 

at higher frequencies. The distribution of the rest-frame synchrotron peak frequency (Vp^^j^) in the spectral energy distribution (SED) of FSRQs 

is the same in all the blazar samples with ( v^^^,) = IQ'^ Hz, while the mean inverse Compton peak frequency, ^v'-^,^^^), ranges from 10^' to 
10^^ Hz. The distributions of v^,^^, and v^^^ of BL Lacs are much broader and are shifted to higher energies than those of FSRQs; their shapes 
strongly depend on the selection method. The Compton dominance of blazars, defined as the ratio of the inverse Compton to synchrotron peak 
luminosities, ranges from less than 0.2 to nearly 100, with only FSRQs reaching values larger than about 3. Its distribution is broad and depends 
strongly on the selection method, with y-ray selected blazars peaking at ~ 7 or more, and radio-selected blazars at values close to 1, thus implying 
that the common assumption that the blazar power budget is largely dominated by high-energy emission is a selection effect. A comparison of 
our multi-frequency data with theoretical predictions shows that simple homogeneous SSC models cannot explain the simultaneous SEDs of most 
of the y-ray detected blazars in all samples. The SED of the blazars that were not detected by Fermi-hAT may instead be consistent with SSC 
emission. Our data challenge the correlation between bolometric luminosity and v^^^ predicted by the blazar sequence. 



1. Introduction 

Blazars are jet-dominated extragalactic objects characterized by 
the emission of strongly variable and polarized non-thermal ra- 
diation across the entire electromagnetic spectrum, from radio 
waves to y-rays (e.g., |Urry & Padovani |[r995 ). As the extreme 
properties of these sources are due to the relativistic amplifica- 
tion of radiation emitted along a jet pointing very close to the 
line o f sight (e.g., [Blandford & Rees|[T978l |Urry & Padovani| 
1995^, they are rare compared to both objects pointing their jets 
at random angles and radio quiet QSOs where the emitted radia- 
tion is due to thermal or reflection mechanisms ultimately pow- 
ered by accretion onto a supermassive black hole (e.g., Abdo 



et al. 2010a I. Despite that, the strong emission of blazars at 



all wavelengths makes them the dominant type of extragalac- 
tic sources in the radio, microwave, y-ray, and TeV bands where 
accretion and other thermal emission processes do not produce 
significant amounts of radiation ( jToffblatti et al.|1998[ |Giommi| 



|& Colafranc esco|2004t[Hartman e t al.ll9 991|Abdo et al.l2010a 
Costamante & Ghisemni|2002[ " Colafrance sco & Gio mmi 2006 
|WeekeS||2008| . For these reasons, blazars are hard to distin- 
guish from other sources at optical and X-ray frequencies, while 
they dominate the microwave and y-ray sky at high Galactic 
latitudes. The advent of the Fermi |A twood et al.||2009i ) and 
P/a»cfcn' |Tauber et al.|2010||Planck C ollaborati on|201 la| l satel- 
lites, which are surveying these two observing windows, com- 
bined with the versatility of the Swift observatory (Gehrels et ah] 
2004a), is giving us the unprecedented opportunity to collect 
multi-frequency data for very large samples of blazars and open 



' Planck (http://www.esa.int/Planck) is a project of the European 
Space Agency - ES A - with instruments provided by two scientific con- 
sortia funded by ESA member states (in particular the lead countries: 
France and Italy) with contributions from NASA (USA), and telescope 
reflectors provided in a collaboration between ESA and a scientific con- 
sortium led and funded by Denmark. 
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the way to a potentially much deeper understanding of the 
physics and demographics of these still puzzling objects. 

Blazars can be categorized by their optical properties and the 
shape of their broad-band spectral energy distributions (SEDs). 
Blazar SEDs always show two broad bumps in the log v - 
log vFy space; the lower energy one is usually attributed to syn- 
chrotron radiation while the more energetic one is attributed 
to inverse Compton scattering. Blazars displaying strong and 
broad optical emission lines are usually called flat-spectrum ra- 
dio quasars (FSRQs), while objects with no broad emission lines 
(i.e., rest-frame equivalent width, EW, < 5A) are called BL 



Lac objects. Padovani & Giommi ( 1995) 1 introduced the terms 
LBL and HBL to distinguish between BL Lacs with low and 
high v alues of the peak f requency of the synchrotron bump 



Abdo et al. 



(2010a I extended this definition to all types 



(Vak)- 

of blazars and defined the terms LSP, ISP, and HSP (correspond- 
ing to low, intermediate, and high synchrotron peaked blazars) 
for the cases where v^^^|^< 10'"^ Hz, lO'-^Hz < v^^^^< 10'^ Hz, 



and Vp^^j^> 10'^ Hz, respectively. In the rest of this paper, we use 
the LSP/ISP/HSP nomenclature. 

It is widely recognized that one of the most effective ways 
of studying the physics of blazars is through the use of multi- 
frequency data that is ideally simultaneous. There are several 



examples of studies following this approach (e.g.. 


Giommi et al. 
1996, Fossati 


1995 


von Montigny et al.|i995, Sambruna et al. 


et al. 


1998| Giommi et al.|2002 [Nieppola et al.|2006 Padovani 


etal. 


2006|l, but in most cases the samples are heterogeneous and 



the data are sparse and non-simultaneous. 

The compilation of simultaneous and well-sampled SEDs re- 
quires the organization of complex multifrequency observation 
campaigns, involving the coordination of observations from sev- 
eral observatories. Such large efforts have been carried out only 
rarely, and almost exclusively on the occasion of large flaring 
events in a few bright and well-known blazars, e.g., 3C 454.3 
( [Giommi et al . '2006', 'Abdo et al. '2009 a] [Vercellone et al.|2009| l, 
Mkn421, (|Donnarumma et al. 2009 [ | Abdo et al.||2011| l, and 
PKS 2155-304 QAharonian et al.|2009l ). 

Significant progress has been made with the publication of a 
compilation of quasi-simultaneous SEDs of a large sample of y- 
ray bright blazars (Abdo et al. '2010a'). This is an important step 
forward from previous compilations, as the sample presented is 
statistically representative of the population of bright y-ray se- 
lected blazars, and the data were quasi-simultaneous, that is col- 
lected within three months of the Fermi-hPC[ observations. 

With Planck, Swift, and Fermi-hPsI simultaneously in or- 
bit, complemented by other space and ground-based observato- 
ries, it is now possible to assemble high-quality multi -frequency 
datasets that allow us to build simultaneous and well-sampled 
broad-band spectra of large and statistically well-defined sam- 
ples of active galactic nuclei (AGNs). 

In this and a companion paper (Planck Collaboration 20Tle|, 
we present the first results of a large cooperative program be- 
tween the Planck, Fermi-hP^I, and Swift satellites and a num- 
ber of ground-based observatories, carried out to collect multi- 
frequency data on large samples of blazars selected using differ- 
ent criteria and observed when the sources lie in the field of view 
of Planck. 

In this paper, we concentrate on blazars selected in the soft 
X-ray, hard X-ray, and y-ray bands. We present the simultaneous 
data, test for flux correlations, and estimate some key parame- 
ters characterizing the SEDs. We then compare the results ob- 
tained for the different samples. Detailed fits to models, variabil- 



ity studies, and more complete theoretical interpretations will be 
presented elsewhere. 

Throughout this paper, we define the radio-to-submillimetre 
spectral index ahy S (v) oc v", and we adopt a ACDM cosmol- 
ogy wifli Hq ^ 70km s ' Mpc', fi„, = 0.27, and = 0.73 
( |Komatsuetal.|2011| l. 



2. Sample selection 

To explore the blazar parameter space from different viewpoints, 
we used several diff'erent criteria to select the blazars to be ob- 
served simultaneously by Planck, Swift, and Fermi. In this paper, 
we considered three samples of blazars that are flux-limited in 
the high-energy part of the electromagnetic spectrum: soft X-ray 
(0. 1-2.4 keV) sources f rom the ROSAT A ll-Sky Survey Bright 
Source Catalog (IRXS, Voges et al.|1999 hereafter referred to 
as the RASS sample), hard X-ray (15-150 keV) sources from 
the ^wi/f-BAT 54-month source catalog ( Cusumano et al.|2010 



hereafter referred to as the BAT sample), and y-ray sources from 
the Fer mi-hPX 3-month Bright AGN Source List ( [Abdo et aL 



2009b' hereafter referred to as the FermZ-LAT sample). 

These high-energy-selected samples were complemented by 
a radio flux-limited sample of northern sources (hereafter re- 
ferred to as the radio sample), which is presented in the com- 
panion paper (Planck Collab oFation|20 1 Ie| ) . We used these four 
samples, defined in widely diff'erent parts of the electromagnetic 
spectrum, to try to disentangle the intrinsic properties of blazars 
from the heavy selection eff'ects that often afflict blazar samples. 
In total we considered 175 sources. 

We based our classification of different blazar types on the 
Roma-BZCAT catalog (Ma ssaro et al.|[20T0| , which is a com- 
pilation of known blazars that were carefully checked to deter- 



mine their blazar type in a uniform and reliable way. Massaro 



et al. (20101 divided blazars into three categories: BZQ/FSRQ, 
in which the optical spectrum has broad emission lines; BZB/BL 
Lac objects, in which the optical spectrum is featureless or con- 
tains only absorption lines from the host galaxy; BZU/uncertain 
type, comprising objects for which the authors could not find 
sufficient data to safely determine the source classification, and 



objects that have peculiar characteristics (see Massaro et al. 



2010, for details). According to this classification, 96 of our ob- 
jects are of the FSRQ type, 40 are BL Lacs, and the rest are 
of uncertain type. About 160 were observed by Swift simulta- 
neously with Planck, mostly by means of dedicated target of 
opportunity (ToO) pointings. In the following we describe the 
selection criteria for each high-energy selected sample. Details 
of the radio-selected sample are given in [Planck Collaboration! 
pOlle^ . 

2. 1 . The issue of blazar classification 

The classification of blazars as either featureless (BL Lacs) or 
broad-lined objects (FSRQs), although very simple in principle, 
is neither unambiguous, nor robust. The borderUne between the 
two blazar subclasses, namely 5 A in the source rest-frame for 
the EW of any emission line, was originally defined as a result 
of the optical identification campaigns of the sources discovered 
in the first well-defined and complete samples of (bright) ra- 
dio and X-ray selected objects ( Stickel et al.|199l" Stocke et al 



1991). However, we now know that well-known BL Lac objects 
such as OJ 287 - and BL Lac itself - exhibit emission lines with 
EWs well above the 5 A limit on some occasions (Vermeulen] 
|erari[T995| [Corbett et aL][T996| . Several other BL Lac objects 
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have strong emission lines with EWs just below, and sometimes 
above the 5 A threshold, depending on the variable continuum 
level (see e.g. [Lawrence et al. 1996jiGhisellini et al.|2011, and 
references therein). Well-known FSRQs such as 3C 279 also ap- 
pear nearly featureless during bright states (Plan et al. 1999 1. 
The detection of broad Lyman-o' emission in the UV spectrum 
of classical BL Lacs such as Mkn421 andMknSOl (Stocke et al.l 



|201 1) contributes to the blurring of the distinction between the 
two types of blazars. 

It is difficult to differentiate between BL Lac objects and ra- 
dio galaxies as BL Lacs have been defined as sources for which 
the 4000 A Ca H&K break (a stellar absorption feature in the 
host galaxy) is diluted by non-thermal radiation more than a 
certain amount that was first quantified by Stocke et al.| ( 1991| l 



T996jl, 



Landt et al. (2002i, 



and then revised by Marcha et al. 

and Landt et al. ( 2004[ l. The level of non-thermal blazar light 
around 4000 A reflects the intrinsic radio power of the jet; it can 
be highly variable and depends strongly on the position of the 
peak of the synchrotron emission, thus ensuring that the border 
between BL Lacs and radio galaxies is quite uncertain. |Giommi| 
et al. ( 2012| l tackled the problem of blazar classification using 
extensive Monte Carlo simulations and showed that the observed 
dififerences can be interpreted within a simple scenario where 
FSRQs and BL Lacs share the basic non-thermal emission prop- 
erties. 

In the present study, we relied on the blazar classification 
given in the Roma-BZCAT catalog ( Massaro et al.|20l"0 i, which 
re-assessed the blazar subclass of each object after a critical re- 
view of the optical data available in the literature and in large 
public databases such as the SDSS ( [York et aL]|2000 l. Despite 
that, some uncertainties remain, which may in turn influence our 
conclusions about the differences between BL Lacs and FSRQs. 
However, the large size of our samples ensures that a few mis- 
classifications should not significantly affect our results. To as- 
sess the impact of both blazar misclassification and transitional 
objects in a quantitative way, it is necessary to perform detailed 
simulations. 



2.2. The Fermi-LAT (j-ray flux-limited) sample 

Our y-ray flux-limited sample was created from the Fermi LAT 
Bright Source Lisf]( Ab do et al. 2009 b). We selected afl the high 
Galactic latitude {\b\ > 10°) blazars detected with high signif- 
icance (TS > lOOjj To reduce the size of the sample and en- 
sure that all the sources are well above the Planck sensitivity 
limit for one survey, we considered only the sources with ra- 
dio flux density (taken from BZCAT) S 5 ghz > 1 Jy- We real- 
ized that this is a double cut, with a TS limit at y-ray energies 
and a flux-density limit in the radio band. A TS limit translates 
into different y-ray flux limits depending on the y-ray spec- 
tral slope, with higher sensitivity to flat-spectrum sources (see 
fig. 7 of Abdo et al.„2009b ). Hence, for our statistical consid- 
erations we also considered the subs ample with a flux cut of 
F(E > 100 MeV) > 8 x 10"'^ phcm^^s"', which removed this 
dependence on the spectral slope. 

The sample so defined includes 50 sources, 40 of which 
are brighter than the y-ray flux limit of 8 x 10"** phcm"^ s"'. 
Relaxing the radio flux density cut would have provided a purely 
y-ray flux-limited sample and increased the number of sources 



^ http : //www . asdc . asi . it/f ermibsl/ 

^ The test statistic (TS) is defined as TS = -21n(Lo/Li) with Lq the 
likelihood of the null-hypothesis model and Lj the likelihood of a com- 
petitive model (see e.g. jAbdo et al.j20l0c| . 



to a; 70, but with about 40-50% of the objects with 5 5 ghz < 1 Jy 
being undetected by Planck. 

Details are presented in Table [T] where column 1 gives the 
source common name, column 2 gives the Fermi-LAI name as 
it appears in 'Abdo et al. ( 2010b) l, columns 3 and 4 give the 
source position in equatorial coordinates, columns 5, 6, and 7 
give the redshift, V magnitude, and X-ray flux (0. 1-2.4 keV) 



from BZCAT ( [Massaro et al.|2010| ), col umn 8 gives the 1.4 G Hz 
or 843 MHz fl ux density from NV SS ( jCondon et aL[[T998| or 
from SUMSS ( |Mauch et ani2003]l with Dec < -40°, column 



Abdo et al 



(2010bb r and column 



9 gives the y-ray flux from 

10 gives the date of the Swift ToO observation made when the 
source was within the Planck field of view. 



2.3. The Swift/BAT (hard X-ray flux-limited) sample 

We defined our hard X-ray flux-limited sample starting from th e 
Swift-BAI 54 month source catalog (Cusumano et al. 2010 1, 
and selecting all the sources identified with blazars with X-ray 
flux > 10"" erg cm"^ s"' in the 15-150 keV energy band. The 
BAT catalog includes 70 known blazars that satisfy the X-ray 
flux condition, but many of them are too faint to be detected at 
millimetre wavelengths by Planck. Therefore, although a pure 
X-ray selection would be preferable, we have decided to add 
a mild radio flux-density constraint (^sghz > 100 mJy, with 
S 5 GHz taken from BZCAT) to select only sources that can be 
detected by Planck or for which Planck will be able to provide 
meaningful upper limits, leaving enough sources to build a sta- 
tistically sizable sample. The list of the 34 sources included in 
this sample is given in Table |2] The column description is the 
same as for Table [T] 



2.4. The ROSAT/RASS (soft X-ray flux-limited) sample 

The soft X-ray flux-limited sample was defined starting from the 
RASS catalog (IRXS) ( |Voges et al.[T999| , and selecting aU the 
sources identified with blazars with count rates higher than 0.3 
counts/s in the 0. 1-2.4 keV energy band, and radio flux densities 
(taken from BZCAT) of 5 5 ghz > 200 mJy. The reasons for using 
an additional radio flux constraint are the same as for the hard 
X-ray flux-limited sample, where, however, we chose 200 mJy 
to reduce the size of the sample to be comparable to those of the 
y-ray and hard X-ray samples. We realize that this is a stringent 
cut that removes about two thirds of the sources from the purely 
soft X-ray selected sample. However, all the sources below the 
radio threshold are HSP BL Lacs, thus implying that the sub- 
sample of LBL sources remains purely X-ray flux-limited and, 
consequently, that high v^^^ objects are strongly underrepre- 
sented. The list of the 43 sources included in this sample is given 
in Table |3] The column description is the same as for Tables [T] 
and|2l 



2.5. The radio flux-density limited sample 

The radio flux-density limited sample is presented in the com- 
panion Planck paper ( [Planck Collaboratio n 201 le), where all 
the observational details are given. The sample consists of 104 
bright northern and equatorial radio-loud AGN characterized by 
Sji GHz > 1 Jy as measured with the Metsahovi radio telescope. 



We give average y-ray fluxes from the 1-year Fermi catalog rather 
than the three-month fluxes that were used to define the sample for con- 
sistency with Tables[2]and[3l 

5 



jhttp : //ww w . asdc . asi . it/bat54/ 



3 



p. Giommi et al.: Simultaneous Planck, Swift, and Fermi observations of blazars 
Table 1. The Fermi-hP^ (y-ray TS/flux-limited) sample. 
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Source name 


Fermi-\^M name (IFGL) 


(Jzuuu; 


1 10AAA\ 
(JZUUUJ 


z 


Rmag 


b 


mjy 




Swift obs. date 


Blazar type 


lJy0118-272 


IFGLJ 0120.5-2700 


01 20 31.6 


-27 01 24 


0.557 


15.5 


0.72 


934 


3.7 ± 0.4 




BLLac - LSP 




IFGLJ 01 37.0+4751 


01 36 58.5 


47 5 1 29 


0.859 


17.6 


1.04 


1138 


9.6 ± 0.6 


2010-02-05 


FSRQ -LSP 


PKS 0202-17 


1 FGLJ 0205.0-1702 


02 04 57.6 


-17 01 18 


1.740 


17.3 


0.57 


1220 


1.5 ± 0.3 


2010-01-08 


FSRQ - LSP 


PKS 0208-512 


IFGLJ 021 0.6-5 101 


02 10 46.2 


-51 01 01 


1.003 


14.8 


0.75 


3493 


7.1 ± 0.6 


2009-1 1-26 


Uncertain - LSP 


PKS 0215+015 


IFGLJ 0217.9+0144 


02 17 48.9 


01 44 49 


1.715 


18.7 


2.56 


751 


6.0 ± 0.5 




FSRQ - LSP 


IJy 0218+357 


IFGLJ 0221.0+3555 


02 21 05.5 


35 56 14 


0.944 


20.0 


0.85 


1707 


6.4 ± 0.5 


2010-08-19 


Uncertain - LSP 


4C 28.07 


IFGLJ 0237.9+2848 


02 37 52.4 


28 48 09 


1.213 


18.8 


0.58 


2197 


3.7 ± 0.4 


2010-02-05 


FSRQ - LSP 


PKS 0235+164 


IFGLJ 0238.6+ 1637 


02 38 38.9 


16 36 59 


0.940 


18.5 


1.24 


1941 


32.7 ±1.1 


2010-01-30 


BLLac - LSP 


PKS 0332-403 


IFGLJ 0334.2-4010 


03 34 13.6 


-40 08 25 




17.5 


0.73 


1042 


3.8 ± 0.4 


2010-01-17 


BL Lac - LSP 


PKS 0420-01 


IFGLJ 0423.2-01 18 


04 23 15.8 


-01 20 33 


0.916 


16.7 


1.39 


2726 


5.6 ± 0.5 


2009-08-27 


FSRQ - LSP 


PKS 0426-380 


IFGLJ 0428.6-3756 


04 28 40.4 


-37 56 19 


1 .030 


16.3 


0.42 


753 


25.7 ± 1.0 


2010-08-17 


BLLac - LSP 


PKS 0454-234 


IFGLJ 0457.0-2325 


04 57 03.1 


-23 24 52 


1 .003 


17.9 


0.91 


1727 


32.5 ± 1.1 


2010-02-25 


FSRQ - LSP 


PKS 0528+134 


IFGLJ 0531.0+1331 


05 30 56.4 


13 31 55 


2.070 


18.9 


0.80 


1556 


4.0 ± 0.5 


2009-09-24 


FSRQ - LSP 


PKS 0537-441 


IFGLJ 0538.8-4404 


05 38 50.3 


-44 05 08 


0.892 


16.0 


2.10 


3729 


21.3 ± 0.9 


2010-03-03 


BL Lac - LSP 


PKS 0735+17 


IFGLJ 0738.2+ 1741 


07 38 07.3 


17 42 19 


0.424 


14.5 


0.97 


2258 


4.4 ± 0.5 


2010-10-07 


BL Lac - ISP 


S4 08 14+425 


IFGLJ 0818.2+4222 


08 18 16.0 


42 22 45 


0.530 


19.6 


0.32 


1091 


8.7 ± 0.6 


2010-10-15 


BLLac - LSP 


OJ535 


IFGLJ 0825.0+5555 


08 24 47.2 


55 52 42 


1.417 


18.1 


0.64 


1449 


0.9 ±0.3 


2010-03-28 


FSRQ - LSP 


PKS 085 1+202 


IFGLJ 0854.8+2006 


08 54 48.8 


20 06 30 


0.306 


14.4 


1.72 


1512 


2.7 ± 0.4 


2010-04-10 


BLLac - LSP 


S4 09 17 +44 


IFGLJ 0920.9+4441 


09 20 58.4 


44 41 54 


2.190 


19.2 


1.04 


1017 


14.0 ± 0.7 


2009-10-29 


FSRQ - LSP 


4C55.17 


IFGLJ 0957.7+5523 


09 57 38.1 


55 22 57 


0.896 


16.8 


0.51 


3079 


10.5 ± 0.6 


2009-1 1-01 


FSRQ - LSP 


4C 01.28 


IFGLJ 1058.4+0134 


10 58 29.6 


01 33 58 


0.888 


17.6 


1.08 


3220 


7.1 ± 0.6 


2009-12-03 


Uncertain - LSP 


PKS 1057-79 


IFGLJ 1058.1-8006 


10 58 43.3 


-80 03 54 


0.581 


17.3 


0.43 


534 


2.2 ± 0.4 


2010-08-30 


BL Lac -LSP 


PKS 1 127-145 


IFGLJ 1130.2-1447 


11 30 07.0 


-14 49 27 


1.184 


16.0 


1.39 


5622 


2.4 ± 0.4 


2009-12-28 


FSRQ - LSP 


PKS 1 144-379 


IFGLJ 1 146.9-3812 


11 f\1 O 

11 47 01.3 


-38 12 11 


1.048 


15.7 


0.91 


1804 


2.4 ± 0.4 


2010-06-24 


BL Lac - LSP 


4C 29.45 


IFGLJ 1 159.4+2914 


11 59 31.8 


29 14 44 


0.729 


16.4 


0.84 


2031 


5.3 ± 0.5 


2010-05-28 


FSRQ - LSP 


ON 231 


IFGLJ 1221.5+2814 


12 21 31.6 


28 13 58 


0.102 


14.3 


1.30 


732 


6.9 ± 0.5 


2009-12-10 


BLLac - ISP 


3C273 


1 FGLJ 1 229 . 1 +0203 


12 29 06.7 


02 03 08 


0.158 


14.1 


63.11 


54991 


9.6 ± 0.6 


2010-01-11 


FSRQ - LSP 


PKS 1244-255 


IrOLJ 1246.7-2545 


1 A £. A ^ O 

12 46 46.5 


'^C ACt 

-25 47 49 


0.635 


16.7 


1.30 


1 165 


8.1 ± 0.6 


2010-01-25 


FSRQ - LSP 


3C 279 


IrCjLJ 1256.2-U547 


12 56 1 1.1 


-05 47 21 


0.536 


15.0 


20.90 


971 1 


32.4 ±1.1 


2010-01-15 


FSRQ - LSP 


IJy 1308+326 


IFGLJ 1310.6+3222 


13 10 28.6 


32 20 43 


0.997 


19.6 


0.53 


1687 


6.8 ± 0.5 


2009-12-12 


Uncertain - LSP 


PKS 1502+106 


IFGLJ 1504.4+1029 


15 04 24.9 


10 29 39 


1.839 


19.5 


0.16 


1774 


67.0 ± 1.6 


2010-07-29 


FSRQ - LSP 


4C -05.64 


IFGLJ 151 1.1-0545 


15 10 53.5 


-05 43 07 


1.I9I 


16.9 


1.16 


3569 


2.1 ± 0.4 




FSRQ - LSP 


APLib 


IFGLJ 1517.8-2423 


15 17 41.8 


-24 22 19 


0.048 


12.6 


1.05 


2042 


5.6 ± 0.5 


2010-02-20 


BLLac - LSP 


OS -237.8 


IFGLJ 1625.7-2524 


16 25 46.8 


-25 27 38 


0.786 


20.6 




2521 


4.4 ± 0.6 


2010-08-26 


Uncertain - LSP 


4C 38.41 


IFGLJ 1635.0+3808 


16 35 15.4 


38 08 04 


1.814 


17.3 


0.17 


2726 


6.8 ± 0.5 


2010-03-07 


FSRQ - LSP 


NRAO 512 


IFGLJ 1642.5+3947 


1 ^ Art £1 

16 40 29.6 


39 46 46 


1.660 


17.5 


0.15 


976 


5.6 ± 0.5 


2010-08-06 


FSRQ - LSP 


Mkn501 


IFGLJ 1653.9+3945 


16 53 52.2 


39 45 36 


0.033 


8.3 


36.90 


1558 


8.3 ± 0.6 


'^A 1 t\ t\'^ 1 

2010-03-21 


BL Lac - HSP 


OT081 


IFGLJ 1751.5+0937 


17 51 32.8 


09 39 00 


0.322 


17.0 


1.18 


623 


6.4 ± 0.6 


2010-04-01 


Uncertain - LSP 


S5 1803+784 


IFGLJ 1800.4+7827 


18 00 45.6 


78 28 04 


0.680 


14.7 


0.79 


2223 


3.0 ± 0.4 


2009-10-13 


BL Lac - LSP 




1 Pni T 1 CII 1 9 9007 


1 Q 1 1 OQ ^1 


90 Oft 5S 


1119 


18 9 


1 77 


2714 


4 5 ± 5 


9nOQ 1 04 




PMNJ 1923-2104 


IFGLJ 1923.5-2104 


19 23 32.1 


-21 04 33 


0.874 


16.6 


0.77 


3167 


11.9 ± 0.7 


2010-09-30 


FSRQ - LSP 


IJy 2005-489 


IFGLJ 2009.5-4849 


20 09 25.3 


-48 49 53 


0.071 


II.O 


33.24 


1282 


5.0 ± 0.5 


2009-10-05 


BLLac -HSP 


PKS 2052-47 


IFGLJ 2056.3-47 14 


20 56 16.3 


-47 14 47 


1.491 


18.3 


0.56 


2223 


4.6 ± 0.5 


2010-10-18 


FSRQ - LSP 


S32141 + 17 


IFGLJ 2 143 .4+ 1742 


21 43 35.5 


17 43 48 


0.213 


14.4 


0.63 


651 


4.9 ± 0.5 


2009-11-20 


FSRQ - LSP 


IJy 2144+092 


1FGLJ2147.2+0929 


21 47 10.1 


09 29 46 


1.113 


16.9 


1.49 


934 


4.1 ±0.4 


d 


FSRQ - LSP 


BLLac 


IFGLJ 2202.8+4216 


22 02 43.2 


42 16 40 


0.069 


12.5 


1.58 


6051 


7.1 ±0.6 


2009-12-23 


BLLac -LSP 


PKS 2204-54 


IFGLJ 2207.8-5344 


22 07 43.7 


-53 46 33 


1.215 


18.2 


0.52 


1526 


1.6 ±0.3 


2010-05-03 


FSRQ - LSP 


PKS 2227-08 


IFGLJ 2229.7-0832 


22 29 40.0 


-08 32 54 


1.560 


16.8 


8.74 


968 


4.6 ± 0.5 


2009-11-19 


FSRQ - LSP 


4C 11.69 


IFGLJ 2232.5+ 1144 


22 32 36.4 


1 1 43 50 


1.037 


16.5 


1.26 


7202 


4.1 ±0.4 


2009-11-29 


FSRQ - LSP 


3C 454.3 


IFGLJ 2253.9+ 1608 


22 53 57.7 


16 08 53 


0.859 


13.2 


7.80 


12657 


46.2 ± 1.3 


2009-12-14 


FSRQ - LSP 



Notes. 843 MHz flux is reported for sources with Dec < -40° ; Units of 10 erg cm ^ s ' ; Units of 10 ph cm ^ s ' 
Swift simultaneous observation not available. Optical spectrum completely featureless or not available, redshift unknown 



Although the samples are defined by different criteria, four 
sources are common to all samples. These are the well-known 
objects 3C273, 3C279, MknSOl, and 3C 454.3, which are 
among the brightest objects across the entire electromagnetic 
spectrum. A summary of the number of sources common to more 
than one sample is given in Table]?] 



3. Data analysis 

3. 1 . Ground-based follow-up data 

Following the launch of Planck, several follow-up programs 
with ground-based facilities started collecting simultaneous ra- 
dio and optical data. In this paper, we used data from the obser- 
vatories listed in Table |5] 



3.1.1. APEX 

Some sources from our sample were observed in the submil- 
limetre domain with the 12-m Atacama Pathfinder Experiment 
(APEX) in Chile. The observations were made using the 
LABOCA bolometer array centered at 345 GHz. Data were 
taken at two epochs: September 3^, 2009, and November 12, 
2009. The data were reduced using the script package mini- 
crusl^ version 30-Oct-2009. Uranus was used as a calibrator 
of the flux densities. 

3.1.2. ATCA-PACO 

The Planck-ATCA Co-eval Observations (P ACO) project 
dMassardi et~aL]|2011a| [Bonavera et al.||2011| ) observed 480 
sources selected from the Australia Telescope 20 GHz catalogue 

' http : //www . subirnn . caltech . edu/~ share/ crush/ 
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Table 2. The 5w!/f-BAT (hard X-ray flux-limited) sample. 



X-ray Flux Flux Density F(?rm/-LAT Flux 







R.A. 


Dec. 






0. 1-2.4 keV 


1.4 GHz" 


1-100 GeV 






Source name 


F(?rm/-LAT name 


(J2000) 


(J2000) 


z. 


Rmag 




mJy 




Swift obs. date 


Blazar type 


IIIZW2 




00 10 31.0 


10 58 29 


0.089 


13.9 


6.14 


98 




2010-07-08 


FSRQ - LSP 


S50014+813 




00 17 08.4 


81 35 08 


3.387 


15.9 


0.77 


693 




2010-09-21 


FSRQ - ISP 


lES 0033+595 


IFGLJ 0035. 9+5951 


00 35 52.6 


59 50 03 




17.2 


5.41 


147 


3.2 + 0.5 


d 


BLLac -HSP 


Mkn348 




00 48 47.1 


31 57 25 


0.015 


9.3 


0.12 


292 




ci 


BLLac -LSP 


IJy 0212+735 


IFGLJ 0217.8+7353 


02 17 30.8 


73 49 32 


2.367 


18.8 


0.54 


2272 


1.0 + 0.4 


2010-09-11 


Uncertain - LSP 


NGC 1275 


IFGLJ 0319.7+4130 


03 19 48.1 


41 30 42 


0.018 


12.3 


197.92 


22830 


17.3 + 0.8 


2010-08-09 


FSRQ -LSP 


NRAO 140 


IFGLJ 0334.2+3233 


03 36 30.1 


32 18 29 


1.259 


16.6 


2.80 


2677 


1.0 + 0.4 


2010-08-25 


Uncertain - LSP 


3C120 




04 33 11.0 


05 21 15 


0.033 


13.8 


22.68 


3440 




2010-02-25 


Uncertain - LSP 


PKS 0521-36 


IFGLJ 0522.8-3632 


05 22 57.9 


-36 27 30 


0.055 


11.6 


10.42 


15620 


2.9 ± 0.4 


2010-03-05 


Uncertain - LSP 


PKS 0528+134 


IFGLJ 0531.0+1331 


05 30 56.4 


13 31 55 


2.070 


18.9 


0.80 


1556 


4.0 ± 0.5 


2009-09-24 


FSRQ - LSP 


IJy 0537-286 


IFGLJ 0539. 1-2847 


05 39 54.2 


-28 39 55 


3.104 


19.0 


0.83 


862 


0.9 ± 0.0 


2010-03-12 


FSRQ - LSP 


PKS 0548-322 




05 50 40.6 


-32 16 20 


0.069 


13.1 


26.34 


344 




2010-03-12 


BLLac -HSP 


B2.2 0743+25 


IFGLJ 0746.6+2548 


07 46 25.8 


25 49 02 


2.979 


19.2 


0.38 


417 


0.7 ± 0.2 


2010-10-15 


FSRQ - LSP 


4C 71.07 


IFGLJ 0842.2+7054 


08 41 24.3 


70 53 42 


2.218 


16.8 


5.52 


3823 


1.2 + 0.3 


2010-03-21 


FSRQ - LSP 


Mkn421 


IFGLJ 1104.4+3812 


11 04 27.3 


38 12 31 


0.030 


8.3 


180.94 


767 


26.1 ± 1.0 


2009-11-17 


BLLac -HSP 


PKS 1127-145 


IFGLJ 1130.2-1447 


11 30 07.0 


-14 49 27 


1.184 


16.0 


1.39 


5622 


2.4 ± 0.4 


2009-12-28 


FSRQ - LSP 


PKS 1219+04 


IFGLJ 1222.5+0415 


12 22 22.5 


04 13 15 


0.967 


17.1 


0.91 


800 


0.9 ± 0.3 


2010-07-17 


FSRQ - LSP 


3C273 


IFGLJ 1229.1+0203 


12 29 06.7 


02 03 08 


0.158 


14.1 


63.11 


54991 


9.6 ± 0.6 


2010-01-11 


FSRQ - LSP 


3C279 


IFGLJ 1256.2-0547 


12 56 11.1 


-05 47 21 


0.536 


15.0 


20.85 


9711 


32.4+ 1.1 


2010-01-15 


FSRQ - LSP 


APLib 


IFGLJ 1517.8-2423 


15 1741.8 


-24 22 19 


0.048 


12.6 


1.05 


2042 


5.7 + 0.5 


2010-02-20 


BLLac -LSP 


Mkn501 


IFGLJ 1653.9+3945 


16 53 52.2 


39 45 36 


0.033 


8.3 


36.93 


1558 


8.3 ± 0.6 


2010-03-21 


BLLac -HSP 


ARP102B 




17 19 14.4 


48 58 49 


0.024 


9.4 


0.59 


145 




2010-03-31 


Uncertain - ISP 


PKSB 1830-210 


IFGLJ 1833.6-2103 


18 33 39.8 


-21 03 39 


2.507 


16.6 


0.69 


10896 


10.7 ± 0.8 


2010-09-23 


FSRQ - LSP 


OV -236 


IFGLJ 1925.2-2919 


19 24 51.0 


-29 14 30 


0.352 


17.3 


2.42 


13387 


1.4 + 0.4 


2010-09-30 


FSRQ - LSP 


lES 1959+650 


IFGLJ 2000.0+6508 


19 59 59.8 


65 08 54 


0.047 


11.9 


35.28 


250 


6.0 ± 0.5 


2009-09-26 


BLLac -HSP 


IJy 2126-158 




21 29 12.1 


-15 3841 


3.268 


16.5 


1.54 


590 




2010-05-03 


FSRQ - ISP 


4C 06.69 


IFGLJ 2148.5+0654 


21 48 05.4 


06 57 38 


0.999 


15.1 


1.46 


2589 


0.7 ± 0.3 


2009-11-21 


FSRQ - LSP 


PKS 2149-307 




21 51 55.5 


-30 27 53 


2.345 


17.4 


4.80 


1243 




2010-05-11 


FSRQ - LSP 


BLLac 


IFGLJ 2202.8+4216 


22 02 43.2 


42 16 40 


0.069 


12.5 


1.57 


6051 


7.1 ±0.6 


2009-12-23 


BLLac -LSP 


4C 31.63 




22 03 14.9 


31 45 38 


0.295 


14.3 


3.22 


2878 




2009-11-27 


FSRQ - LSP 


NGC 7213 




22 09 16.2 


-47 10 00 


0.006 


10.3 


35.34 


98 




2010-10-23 


Uncertain 


4C 11.69 


IFGLJ 2232.5+1144 


22 32 36.4 


1 1 43 50 


1.037 


16.5 


1.26 


7202 


4.1 ±0.4 


2009-11-29 


FSRQ - LSP 


3C 454.3 


IFGLJ 2253.9+ 1608 


22 53 57.7 


16 08 53 


0.859 


13.2 


7.80 


12657 


46.2+ 1.3 


2009-12-14 


FSRQ - LSP 


PKS 2325+093 


IFGLJ 2327.7+0943 


23 27 33.5 


09 40 09 


1.843 


18.8 


0.73 


741 


3.0 ± 0.4 


2010-07-18 


FSRQ - LSP 



Notes. 843 MHz flux is reported for sources wi til Dec < -40°. Units of 10 '-ergcm ^s ' ; Units of 10 '"phcm "s ' 
Swift simultaneous observation not available. Optical spectrum completely featureless or not available, redshift unknown 



(AT20G, [Massardi et al.||2011b| l, with flie Australia Telescope 
Compact Array (ATCA) in the frequency range 4.5^0 GHz, at 
several epochs close in time to the Planck observations in the pe- 
riod July 2009 to August 2010. The PACO sample is a complete, 
flux-density limited, and spectrally selected sample of southern 
sources, with the exception of the region with Galactic latitude 
\b\ < 5°. A total of 147 PACO point-like sources have at least 
one observation within ten days of the Planck observations. 

3.1.3. Effelsberg and IRAM 

Quasi-simultaneous cm/mm radio spectra for a larger number of 
Planck blazars were obtained within the framework of a Fermi 



monitoring program of y-ray blazars (F-GAMMA: Fuhrmann 
[eTaL] [20571 [Angelakis et"aL][2008| on the Effelsberg 100-m 
and IRAM 30-m telescopes. The frequency range was 2.64- 
142 GHz. 

The Effelsberg measurements were conducted with the sec- 
ondary focus heterodyne receivers at 2.64, 4.85, 8.35, 10.45, 
14.60, 23.05, 32.00, and 43.00 GHz. The observations were per- 
formed quasi-simultaneously with cross-scans, that is by slew- 
ing over the source position in azimuth and elevation with the 
number of sub-scans chosen to reach the desired sensitivity 
(for details, see |Fuhrmann et al.|[2008 Angelakis et al.|[2008[ l. 
Pointing offset, gain, atmospheric opacity, and sensitivity cor- 
rections were applied to the data. 

The IRAM 30-m observations were carried out with cal- 
ibrated cross-scans using the EMIR horizontal and vertical 
polarization receivers operating at 86.2 GHz and 142.3 GHz. 



The opacity-corrected intensities were converted into the stan- 
dard temperature scale and finally corrected for small re- 
maining pointing offsets and systematic gain-elevation effects. 
Conversion to a standard flux density scale was based on fre- 
quent observations of primary calibrators (Mars, Uranus) and 
secondary calibrators (W3(OH), K3-50A, NGC 7027). 

From this program, radio spectra measured quasi- 
simultaneously with the Planck observations were collected for 
a total of 37 Planck blazars during the period August 2009 to 
June 2010. Results are reported in Table |6] and [7] 

3.1.4. Medicina 

The Simultaneous Medicina Planck Experiment (SiMPlE, 
Procopio et al. 2011) used the 32-m Medicina single dish to 



make almost simultaneous observations at 5 GHz and 8.3 GHz 
of the 263 sources of the NEWPS sample ( Massardi et al.|2009| l 
with Dec > 0°, and partially overlapping with the PACO obser- 
vations for -10° < Dec < 0°. The project began in June 2010 
and finished in June 2011, observing our sample several times 
throughout two complete Planck surveys. It does not overlap 
with the Planck first survey. 



3.1.5. Metsahovi 

The 37 GHz observations were made with the 13.7-m Metsahovi 
radio telescope using a 1 GHz bandwidth, dual-beam receiver 
centered at 36.8 GHz. We performed ON-ON observations, by 
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Table 3. The ROSATIRASS (soft X-ray flux-limited) sample. 















X-ray Flux 


Flux Density 


Fertni-h,AT Flux 










R.A. 


Dec. 






y.l— z.4KeV 


1.4 GHz " 


1—1 Ou ueV 






Source name 


termi-LAl name 


(J 2000} 


(J 2000) 




Rmag 


b 


mjy 




Swift obs. date 


Blazar type 


III ZW 2 




00 10 31.0 


10 58 29 


0.089 


13.9 


6.14 


98 




2010-07-08 


FSRQ - LSP 


GB6J 0214+5145 




02 14 17.9 


51 44 52 


0.049 


16.5 


4.58 


294 






BL Lac - HSP 


3C 120 




04 33 1 1 .0 


05 21 15 


0.033 


13.8 


22.68 


3440 




2010-02-25 


Uncertain - LSP 


PKS 0521-36 


IFGLJ 0522.8-3632 


05 22 57.9 


-36 27 30 


0.055 


1 1.6 


10.42 


15620 


2.9 ± 0.4 


2010-03-05 


Uncertain - LSP 


PKS 0548-322 




05 50 40.6 


-32 16 20 


0.069 


13.1 


26.34 


344 




2010-03-12 


BL Lac - HSP 


IRAS-L 06229-643 




06 23 07.6 


-64 36 20 


0.129 


13.7 


5.34 


274 




2010-08-18 


FSRQ - LSP 


4C 7 1 .07 


IFGLJ 0842.2+7054 


08 41 24.3 


70 53 42 


2.218 


16.8 


5.52 


3823 


1.2 ± 0.3 


2010-03-21 


FSRQ - LSP 


B2 0912+29 


IFGLJ 0915.7+2931 


09 15 52.4 


29 33 24 




15.0 


6.25 


342 


2.1 ± 0.3 


2010-10-28 


BLLac -HSP 


PKS 0921-213 




09 23 38.8 


-21 35 47 


0.053 


12.8 


4.80 


268 




2010-05-02 


Uncertain - LSP 


IH 1013+498 


IFGLJ 1015.1+4927 


10 15 04.1 


49 26 00 


0.212 


15.1 


13.23 


378 


8.7 ± 0.6 


2010-04-24 


BL Lac - HSP 


IRXSJ 105837.5+562816 


IrGLJ 1058.6+5628 


10 58 37.7 


56 28 11 


0.143 


14.0 


3.13 


228 


5.7 ± 0.5 


'^n 1 A A /I 1 o 

2010-04-18 


BL Lac - HSP 


PKS 1 124-186 


IFGLJ 1 126.8-1854 


1 1 27 04.3 


-18 57 17 


1.048 


19.2 


5.33 


536 


2.4 ± 0.4 


2010-06-10 


FSRQ - LSP 


B2 1128+31 




11 31 09.4 


31 14 05 


0.289 


15.8 


5.02 


370 




2009-1 1-28 


FSRQ - HSP 


S5 1133+704 


IFGLJ 1 136.6+7009 


1 1 36 26.4 


70 09 27 


0.045 


11.0 


35.08 


327 


1.7 ± 0.3 


2009-10-27 


BL Lac - HSP 


4C 49.22 




11 53 24.4 


49 3 1 08 


0.334 


16.9 


3.31 


1572 




20091 1-17 


FSRQ - LSP 


ON 325 


IFGLJ 1217.7+3007 


12 17 52.0 


30 07 00 


0.130 


14.5 


24.90 


572 


6.7 ± 0.6 


2009-12-03 


BL Lac - HSP 


PKS 1217+02 




12 20 1 1.8 


02 03 42 


0.241 


15.6 


2.78 


672 




2010-06-24 


FSRQ - LSP 


3C 273 


IFGLJ 1229.1+0203 


12 29 06.7 


02 03 08 


0.158 


14.1 


63.11 


54991 


9.6 ± 0.6 


2010-01-11 


FSRQ - LSP 


PG 1246+586 


IFGLJ 1248.2+5820 


12 48 18.7 


58 20 28 




14.5 


3.99 


245 


4.5 ± 0.4 


2010-05-20 


BLLac - ISP 


3C 27y 


IFGLJ 1256.2-0547 


12 56 11.1 


-05 47 2 1 


0.536 


15.0 


20.85 


971 1 


32.4 ±1.1 


2010-01-15 


FSRQ - LSP 


IJy 1302-102 




13 05 33.0 


-10 33 19 


0.286 


14.4 


4.20 


71 1 






FSRQ - ISP 


GB6B 1347+0955 




13 50 22.1 


09 40 10 


0.133 


13.6 


3.74 


293 




2010-07-18 


Uncertain - ISP 


IWGAJ 1407.5-2700 




14 07 29.7 


-27 01 04 


0.022 


9.7 


15.28 


646 




2010-02-12 


Uncertain - HSP 


3C 298.0 




14 19 08.1 


06 28 34 


1.437 


16.4 


0.00 


6100 




2010-07-30 


Radio Galaxy 


BZQJ 1423+5055 




14 23 14.1 


50 55 37 


0.286 


15.1 


3.35 


178 




2010-07-13 


FSRQ - HSP 


PG 1424+240 


IFGLJ 1426.9+2347 


14 27 00.3 


23 48 00 




14.5 


3.57 


430 


10.2 ± 0.6 


2010-01-22 


BLLac - ISP 


IRXSJ 145603.4+504825 




14 56 03.6 


50 48 25 


0.479 


18.1 


13.02 


220 




2009-12-25 


BL Lac - HSP 


BZQJ 1507+0415 




15 07 59.7 


04 15 12 


1.701 


19.0 


6.1 1 


167 




'^Al A AO {\C 

2010-08-05 


FSRQ - LSP 


PG 1553+113 




15 55 43.0 


11 1 1 24 




14.0 


17.85 


312 




2010-02-05 


BLLac - ISP 


WE 1601 + 16W3 




16 03 38.0 


15 54 02 


0.110 


13.4 


4.14 


97 




2010-08-14 


Uncertain - HSP 


3C 345 


IFGLJ 1642.5+3947 


16 42 58.8 


39 48 37 


0.593 


15.0 


2.52 


7099 


5.6 ± 0.5 


2010-03-06 


FSRQ - LSP 


Mkn 501 


IFGLJ 1653.9+3945 


16 53 52.2 


39 45 36 


0.033 


8.3 


36.93 


1558 


8.3 ± 0.6 


2010-03-21 


BL Lac - HSP 


1 PC 1 "7/1 1 Q/^ 


icr^l T yiAA 0_i_lO'i/l 
lrOL,J I /44.Z+iyj4 


in A'l Q 


iy jj uy 


n AC/1 


IZ. / 






1 . 1 ± U.J 


d 


dl i^ac - ijr 


PKS 1833-77 




18 40 38.4 


-77 09 28 


0.018 


8.3 


5.93 


1108 




2010-03-11 


Uncertain - ISP 


lES 1959+650 


IFGLJ 2000.0+6508 


19 59 59.8 


65 08 54 


0.047 


11.9 


35.28 


250 


6.0 ± 0.5 


2009-09-26 


BLLac -HSP 


IJy 2005-489 


IFGLJ 2009.5-4849 


20 09 25.3 


-48 49 53 


0.071 


11.0 


33.24 


1282 


5.0 ± 0.5 


2009-10-05 


BLLac -HSP 


PKS 2149-307 




21 51 55.5 


-30 27 53 


2.345 


17.4 


4.80 


1243 




2010-05-11 


FSRQ - LSP 


NGC7213 




22 09 16.2 


-47 10 00 


0.006 


10.3 


35.34 


98 




2010-10-23 


Uncertain - HSP 


PKS 2227-399 




22 30 40.2 


-39 42 52 


0.318 


16.0 


4.23 


369 




2010-05-09 


Uncertain - ISP 


3C 454.3 


IFGLJ 2253.9+1608 


22 53 57.7 


16 08 53 


0.859 


13.2 


7.80 


12657 


46.2 ± 1.3 


2009-12-14 


FSRQ - LSP 


PKS 2300- 18 




23 03 02.9 


-18 41 25 


0.129 


15.5 


5.16 


861 




2010-05-30 


Uncertain - ISP 


PKS 233 1-240 




23 33 55.2 


-23 43 40 


0.048 


11.4 


3.92 


782 




2010-06-05 


Uncertain - ISP 


lES 2344+514 


IFGLJ 2347. 1 +5 142 


23 47 04.8 


51 42 17 


0.044 


10.7 


7.71 


250 


1.4 + 0.3 


2010-01-17 


BLLac - HSP 



Notes. 843 MHz flux is reported for sources witii Dec < -40°. Units of 10 erg cm ^ s ' ; Units of 10 ph cm - s ' 
Swift simultaneous observation not available. Optical spectrum completely featureless or not available, redshift unknown 



Table 4. Summary of the samples, blazar types, and selection methods considered in this paper. 





Selection 


No. of 


Blazars 






Sources in common 






Sample 


band 


sources 


FS/BL/Unc. 


Other AGN 


RASS 


BAT 


Fermi-hAI 


Radio 


All 


RASS 


Soft X-ray 


43 


15/16/11 


1 




12 


5 


9 


4 


BAT 


Hard X-ray 


34 


21/7/6 




12 




9 


16 


4 


Fermi-LAT * 


y-ray 


50 


28/16/6 




5 


9 




23 


4 


Fermi-LAT FL" 


y-ray 


40 


27/8/5 




3 


7 




19 


3 


Total this paper 




105 


52/32/20 


1 












Radio 


radio 


104 


73/18/10 


3 


9 


16 


23 




4 



• Total Fenni-hAT sample (TS limited), 

** Flux-limited Fermi-LKT sample F(E > 100 Me V) > 8 x lO"'* ph cm"^ s"' 



alternating between the source and the sky in each feed horn. A 
typical integration time for obtaininig one flux density data point 
was 1200-1400 s. The telescope detection limit at 37 GHz was 
~0.2 Jy under optimal conditions. Data points with a signal-to- 
noise ratio (S/N) smaller than four are handled as non-detections. 
The flux-density scale was set by observations of DR 2 1 . Sources 
NGC7027, 3C274, and 3C84 were used as secondary calibra- 
tors. A detailed description of the data reduction and analysis is 



given in ( |Terasranta et al.|["l998| l. The error estimate in the flux 
density includes the contribution from the measurement rms and 
the uncertainty in the absolute calibration. 

3.1.6. OVRO 

Some of the sources in our samples were monitored at 15 GHz 
using the 40-metre telescope of the Owens Valley Radio 
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Table 5. Optical and radio observatories participating in the Planck multi-frequency campaigns. 



Radio observatory 


Frequencies [GHz] 


APEX, Chile 


345 


ATCA, Australia 


4.7, 5.2, 5.8, 6.3, 8.2, 8.7, 9.3, 9.8, 17.2, 17.7, 18.3, 18.8 

23.2, 23.7, 24.3, 24.8, 32.2, 32.7, 33.3, 33.8, 38.2, 38.7, 39.3, 39.8 


Effelsberg, Germany 
IRAM, Spain 
Medicina, Italy 


Z.O^, ^.OJ, o.JJ, 1U.4J, 14. DU, Aj.yjJ, jZ.UU, ^D.\J\l 

86.2, 142.3 

Z), o 


Metsahovi, Finland 


37 


OVRO, USA 


15 


RATAN, Russia 


1.1,2.3,4.8, 7.7, 11.2,21.7 


UMRAO, USA 


4.8, 8.0, 14.5 


VLA, USA 


5, 8, 22, 43 


Optical observatory 


Band 


KVA, Spain 
Xinglong, China 


R 
I 



Observatory as part of a larger monitoring program (Richards 
et al.|20lT l. The flux density of each source was measured ap- 
proximately twice weekly, with occasional gaps due to poor 
weather or instrumental problems. The telescope was equipped 
with a cooled receiver installed at prime focus, with two 
symmetric off-axis corrugated horn feeds that are sensitive 
to left circular polarization. The telescope and receiver com- 
bination produces a pair of approximately Gaussian beams 
(157 arcsec FWHM), separated in azimuth by 12.95 arcmin. The 
receiver has a central frequency of 15.0 GHz, a 3.0 GHz band- 
width, and a noise-equivalent reception bandwidth of 2.5 GHz. 
Measurements were made using a Dicke-switched dual-beam 
system, with a second level of switching in azimuth where we 
alternated between source and sky in each of the two horns. 
Our calibration is refeiTed to 3C 286, for which a flux density 
of 3.44 Jy at 15 GHz is assumed ( |Baars et aL||1977i i. Details of 



the observations, calibration, and analysis are given by ( Richards 
let al.|2011| l. 

3.1.7. RATAN 

A six-frequency broadband radio spectrum was obtained with 
the RATAN-600 radio telescope in transit mode by observing si- 
multaneo usly at 1 .1, 2.3, 4.8, 7.7, 11.2, and 21.7 GHz ( |Parijskij| 
[1993, Berlin & Friedman 1996). Data were reduced using the 
RATAN standard software FADPS (Flexible Astronomical Data 
Processing System) reduction package ( [Verkhodanov 1997[ l. 
The flux density measurement procedure at RATAN-600 is de- 
scribed by ( )Aliakberov et al.,1985) ). 



servations of Cassiopeia A ( |Baars et al. 19771. Details of the 
calibration and analysis techniques are described by Aller et al. 
(1985J . 

3.1.9. VLA 

The Very Large Array (VLA) and, since Spring 2010, the 
Expanded VLA (EVLA), observed a subset of the sources as 
simultaneously as possible. Most of the VLA and EVLA runs 
were performed in one to two hour chunks of time. We ob- 
served during a one-hour chunk of time, in addition to flux cal- 
ibrators and phase calibrators, typically 5-8 Planck sources. In 
many cases, VLA flux-density and phase calibrators were them- 
selves of interest, since they were bright enough to be detected 
by Planck. For these bright sources, the integration times could 
be extremely short. 

Integration times were about 30 s at 4.86 GHz and 8.46 GHz, 
100 s at 22.46 GHz, and 120 s at 43.34 GHz. All VLA/EVLA 
flux density measurements were calibrated using standard val- 
ues for one or both of the primary calibrator sources used by 
NRAO, 3C48 or 3C286, and the u-v data were flagged, cal- 
ibrated, and imaged using standard NRAO software (AIPS or 
CASA). It is important to bear in mind that the VLA and EVLA 
were in different configurations at different times in the several 
months duration of the observations. As a consequence, the an- 
gular resolution changed, becoming, for a given configuration, 
much higher at higher frequencies. For that reason, sources that 
appeared to be resolved in any VLA configuration or at any VLA 
frequency were carefully flagged. 



3.1.8. UMRAO 

Centimetre-band observations were obtained with the University 
of Michigan 26-m prime focus paraboloid equipped with ra- 
diometers operating at central frequencies of 4.8, 8.0, and 
14.5 GHz. Observations at all three frequencies utilized rotating 
polarimeter systems permitting both total flux density and linear 
polarization to be measured. A typical measurement consisted 
of 8 to 16 individual measurements over a 20-40 minute time 
period. Frequent drift scans were made across stronger sources 
to verify the telescope pointing connection curves; and observa- 
tions of program sources were intermixed with observations of 
a grid of calibrator sources to correct for temporal changes in 
the antenna aperture efficiency. The flux scale was based on ob- 



3.1.10. KVA 

Optical observations were made with the 35 cm KVA (Kunliga 
Vetenskapsakademiet) telescope at La Palma, Canary islands. 
All observations were made through the /?-band filter (/leff = 
640 nm) using a Santa Barbara ST-8 CCD camera with a gain 
factor of 2.3 e /ADU and readout noise of 14 electrons. Pixels 
were binned 2x2 pixels giving a plate scale of 0.98 arcsec/pixel. 
We obtained 3-6 exposures of 180 s per target. The images were 
reduced in the standard way of subtracting the bias and dark 
frames and dividing by twilight flat-fields. The fluxes of the tar- 
get and 3-5 stars in the target field were measured with aper- 
ture photometry and the magnitude difference between the target 
and a primary reference star in the same field was determined. 
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The use of differential mode makes the observations insensitive 
to variations in atmospheric transparency and accurate measure- 
ments can be obtained even in partially cloudy conditions. The 
/?-band magnitude of the primary reference star was determined 
from observations made on photometric nights, using compari- 
son stars in known blazar fields as calibrators ( Fiorucci & Tosti| 



1996 ; Fiorucci et al. 1998 ; Raiteri et al.IT998l | Villata et al. 1 1 998 1 
Nilsson et al. 2007 ) and taking into account the color term of the 
iJ-band filter employed. After the /?-band magnitude of the pri- 
mary reference star was determined, the object magnitudes were 
computed from the magnitude differences. At this phase we as- 
sumed V - R - 0.5 for the targets. Several stars in the field were 
used to check the quality of the photometry and stability of the 
primary reference. The uncertainties in the magnitudes include 
the contribution from both measurement and calibration errors. 



3.1.11. Xinglong 

The monitoring at Xinglong Station, National Astronomical 
Observatories of China, was performed with a 60/90 cm f/3 
Schmidt telescope. The telescope is equipped with a 4096x4096 
E2V CCD, which has a pixel size of 12yum and a spatial res- 
olution of r.'3 pixel"' . The observations were made with an 
intermediate-band filter, the / filter Its central wavelength and 
passband width are 6685 A and 5 14 A, respectively. The expo- 
sure times were mostly 120s but ranged from 60s to 180s, de- 
pending on weather and lunar phase. 

3.2. Planck microwave data 

Planck UTauber et al.|2010| [Planck Collaboration||2011a| l is the 
third generation space mission to measure the anisotropy of the 
cosmic microwave background (CMB). It observes the sky in 
nine frequency bands covering 30-857 GHz with high sensitiv- 
ity and angular resolution from 31' to 5'. Full sky coverage is 
attained in about seven months. The Low Frequency Instrument 



(LFI;|Mandolesi et al.|2010l|Bersaneni et al.|2010||Mennella, A.[ 
etal.|2011| l covers the 30, 44, and 70 GHz bands with amplifiers 
cooled to 20 K. The High Frequency Instrument (HFI; |Lamarre| 
let al.|2010| [Planck HFI Core Team|20lTa| l covers the 100, 143, 
217, 353, 545, and 857 GHz bands with bolometers cooled to 
0.1 K. Polarization is measured in all but the highest two bands 
( Leahy et al.|2010 Rosset et al.|2010) . A combination of radia- 
tive cooling and three mechanical coolers produces the temper- 



the second (OD 275-456), and the beginning of the third Planck 
survey (OD 457-550) through the LFI and HFI pipelines de- 
scribed in |Zacchei, A. et aTjpOlll l and [Planck HFI Core Team] 
( [201 lb| l, and we extracted the flux densities from each map us- 
ing IFACMEX, which is an implementation of the Mexican Hat 
Wavelet 2 (MHW2) algorithm available at the LFI DPC. The 
MHW2 tool has been extensively used to detect point-like ob- 
jects in astronomical images, both with simulations from vari- 
ous experiments and data from the WMAP, Planck, and Herschel 



satelU tes ([Gonzalez-Nuevo et al.[J2006[ [Lopez-Caniego et aL 
2006) [2007t [Massardi et al.[|2009{ . This wavelet is defined as 
the fourth derivative of the two-dimensional Gaussian function, 
where the scale of the filter R is optimized to look for the max- 
imum in the S/N of the sources in the filtered map. In practice, 
the IFCAMEX code, our implementation of the MHW algo- 
rithm, deals with flexible image transport system (FITS) maps 
in Healpix format (Go rski et al.|| 2005) and can be used to de- 
tect sources in the whole sky or at the position of known objects. 
For this analysis, we looked for objects above the 4cr level at the 
positions corresponding to the 105 sources of our sample. For 
objects with S/N smaller than four we adopted the 4cr level as 
an upper limit. The results of the analysis of Planck simultane- 
ous data are reported in Table |8] where columns 1 and 2 give the 
source name, columns 3 and 4 give the observation start and end 
times, and columns 5-13 give the flux densities in units of Jy at 
30, 44, 70, 100, 143, 217, 353, 545, and 857 GHz. 

Owing to source variability, we do not expect these simulta- 
neous flux densities to be the same as the time-averaged ERCSC 
measurements, except in the case of the sources that were ob- 
served only once during the ERCSC time range and for which 
we estimated the Planck flux densities, measured simultaneously 
with the Swift observation, in the same period. We verified that, 
for the sources fulfilling these requirements, the flux densities 
extracted for this paper are in good agreement with those of the 
ERCSC. 

In addition to simultaneous Planck data, we also used 
ERCSC flux densities in both our analysis of flux correlations 
(Sect. 7.2 1 and the SED plots described in Sect. [5] 



3.3. Swift optical, UV, and X-ray data 



atures needed for the detectors and optics (Planck Collaboration 



2011b[). Two data processing centers (DPCs) check and calibrate (UVOT, [Roming et al 



The Swift Gamma-Ray-Burst (GRB) Explorer ( [Gehrels et al.[ 
2004b I is a multi-frequency space observatory devoted to the 
discovery and rapid follow-up of GRBs. There are three instru- 
ments on board the spacecraft: the UV and Optical Telescope 
]2005] l, the X-Ray Telescope (XRT, 



the data and make maps of the sky ( [Planck HFI Core Team 
201 lb Zacchei, A. et al.|201 l )l. Planck's sensitivity, angular res- 



olution, and frequency coverage make it a powerful instrument 
for Galactic and extragalactic astrophysics as well as cosmol- 
ogy. Early astrophysics results are given in Planck Collaboration, 
2011h-z. 

The Early Release Compact Source Catalog (ERCSC, 
[Planck Collaboration|201 Ic) contains all sources, both Galactic 
and extragalactic, detected with high confidence over the full 
sky during the period between August 12, 2009 and June 6, 
2010 (corresponding to Planck operational days 91 to 389). 
The ERCSC only contains average intensity information for the 
sources. However, many of the sources were observed more than 
once during the time period spanned by the ERCSC, and some 
of the Swift observations used for this paper were carried out be- 
tween June and October 2010. Therefore, to have simultaneous 
data, we produced independent maps for the first (OD 91-274), 



[Burrows et al.|2005[ ) sensitive to the 0.3-10.0 keV band, and the 
Burst Alert Telescope (BAT, [Barthelmy et al.|2005| l sensitive to 
the 15-150keV band. Although the primary scientific goal of 
the satellite is the observation of GRBs, the wide frequency cov- 
erage is suitable for blazar studies, because it covers the region 
where the transition between synchrotron and inverse Compton 
emission usually occurs. 

When not engaged in GRB observations. Swift is available 
for target of opportunity (ToO) requests, and the Swift team de- 
cided to devote an average of three ToO observations per week 
to this project for simultaneous observations of blazars. 



3.3.1. UVOT 

The Swift UVOT telescope can produce images in each of its 
six lenticular filters (V, B, U, UVWl, UVM2, and UVW2). 
However, in an effort to reduce the use of mechanical parts after 
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several years of orbital operations, observations are caiTied out 
using only one filter, unless specifically requested by the user 
Thus images are not always available for all filters. 

The photometry analysis of all our sources was performed 
using the standard UVOT software distributed within the 
HEAsoft 6.8.0 package and the calibration included in the latest 
release of the "Calibration Database". A specific procedure was 
developed at the ASI Science Data Center (ASDC) to process all 
ToO observations requested for the blazar sample. Counts were 
extracted from apertures of 5" radius for all filters and converted 
to fluxes using the standard zero points ( [Poole et al.|2008| l. The 
fluxes were then de-reddened using the appropriate values of 
E{B - V) for each source taken from (Schlegel et al.,1998j with 
Ai/E{B - V) ratios calculated for UVOT filters using the mean 
Galactic interstellar extinction curve from ("Fitzpatrick 1999 1. 
No variability was detected within single exposures in any filter. 
The processing results were carefully validated including checks 
for possible contamination by nearby objects within the source 
and background apertures. Some sources, such as 3C273 and 
NGC 1275, needed special analysis, and results for some other 
sources had to be discarded. 

The results of the UVOT data analysis are summarized in 
Table |9] where columns 1 and 2 give the source name, columns 
3 and 4 give the observation date and the Swift observation ID, 
and the remaining columns give the magnitudes in the six UVOT 
filters with errors. 

3.3.2. XRT 

The Swift XRT is usually operated in "auto-state" mode, which 
automatically adjusts the CCD read-out mode to the source 
brightness, in an attempt to avoid pile-up ( [Burrows et al.|2005] 



Hill et al. 2004| l. As a consequence, some of the data were 



collected using the most sensitive photon counting (PC) mode, 
while windowed timing (WT) mode was used for bright sources. 

The XRT data were processed with the XRTDAS software 
package (v. 2.5.1, [CapaAi et al.[ ( |2005[ )) developed at the ASI 
Science Data Center (ASDC) and distributed by the NASA High 
Energy Astrophysics Archive Research Center (HEASARC) 
within the HEASoft package (v. 6.9). Event files were calibrated 
and cleaned with standard filtering criteria using the xrtpipeline 
task and the latest calibration files available in the Swift CALDB. 
Events in the energy range 0.3-10 keV with grades 0-12 (PC 
mode) and 0-2 (WT mode) were used for the analysis. 

Events for the spectral analysis were selected within a cir- 
cle of 20 pixels (~ 47") radius, which encloses about 90% of 
the point spread function (PSF) at 1.5 keV ( [Moretti et al.|2005| l, 
centered on the source position. When the source count rate is 
above ~ 0.5 counts/s, the PC mode data are significantly affected 
by pile-up in the inner part of the PSF. In these cases, and after 
comparing the observed PSF profile with the analytical model 
derived by ([Moretti et al.|2005)l, we removed pile-up effects by 
excluding events detected within up to 6 pixels from the source 
position, and used an outer radius of 30 pixels. The value of the 
inner radius was evaluated individually for each observation af- 
fected by pile-up, in a way that depended on the observed source 
count rate. 

Ancillary response files were generated with the xrtmkarf 
task by applying corrections for the PSF losses and CCD de- 
fects. Source spectra were binned to ensure a minimum of 20 
counts per bin when utilizing the minimization fitting tech- 
nique. We fitted the spectra adopting an absorbed power-law 
model with photon index Ff. When deviations from a single 
power-law model were found, we adopted a log-parabolic law 



of the form F(E) = KE^»+b\ogE) ( [Massaro et aL]|2004| l, which 
has been shown to fit the X-ray spectrum of blazars of the HSP 
type well (e.g. Giommi et al. 2005 ; Tramacere et al. 2009). This 
spectral model is described by only two parameters: a, the pho- 
ton index at 1 keV, and b, the curvature of the parabola. For both 
models, the amount of hydrogen-equivalent column density (A^h) 
was fixed to the Galactic value along the line of sight ( [Kalberla] 



etaL 2005). For a fraction of the sources fitted with a power- 



law model (~15%), we found evidence of an absorption excess 
at low energies and the hydrogen column density A^h parameter 
was left free. 

The results of the spectral fits with a power-law model and 
Galactic A^h are shown in Table 10 where columns 1 and 2 give 
the source name, column 3 gives the Swift observation date, col- 
umn 4 gives the best-fit photon index T^, column 5 gives the 
Galactic A^h, columns 6 and 7 give the 0.1-2.4 and 2-10 keV 
X-ray fluxes, column 8 gives the value of the reduced x^i and 
column 9 gives the number of degrees of freedom. 



In Table 1 1 we report data obtained using a log-parabola to 
describe the spectrum model. Columns 1 and 2 give the source 
name, column 3 gives the Swift observation date, columns 4 and 
5 give the log parabola parameters a and b, column 6 gives the 
Galactic A^h^ columns 7 and 8 give the 0.1-2.4 and 2-lOkeV 
X-ray fluxes, column 9 gives the value of the reduced and 
column 10 gives the number of degrees of freedom. 



Finally, in Table 12 we report the results of the spectral fits 
that were performed leaving the hydrogen column density A^h 
to vary as a free parameter. Columns 1 and 2 give the source 
name, column 3 gives the Swift observation date, column 4 gives 
the best-fit photon index F,, column 5 gives the estimated A^h, 
columns 6 and 7 give the 0.1-2.4 and 2-10 keV X-ray fluxes, 
column 8 gives the value of the reduced and column 9 gives 
the number of degrees of freedom. 

3.4. ferm-LAT y -ray data 

The Large Area Telescope (LAT) on-board Fermi is an electron- 
positron pair conversion telescope sensitive to y-rays of ener- 
gies from 20MeV to > 300 GeV The Fermi-LAT consists of 
a high-resolution silicon microstrip tracker, a Csl hodoscopic 
electromagnetic calorimeter, and an anticoincidence detector for 
charged particle background identification. A full description of 
the instrument and its performance can be found in |Atwood et al. 
(2009j ). The large field of view (-2.4 sr) allows the LAT to ob- 
serve the full sky in survey mode every 3hr. The LAT point 
spread function (PSF) depends strongly on both the energy and 
the conversion point in the tracker, but less on the incidence an- 
gle. 

The LAT y-ray spectra of all AGN sources are studied in 
[Abdo et al.[ ( [2010b| l based on 1 1 months of Fermi-LAT data. 
Here we derived the y-ray spectra of the blazars for which we 
built the simultaneous SEDs, integrating for two weeks encom- 
passing the whole duration of the Planck observations. 

The Fen«/-LAT data considered for this analysis cover the 
period from August 4, 2008 to November 4, 2010 and were 
analyzed using the standard FermZ-LAT ScienceTools software 
packag^ (version v9rl6 and selecting for each source only pho- 
tons of energies above 100 MeV belonging to the diffuse class 
(Pass6 V3 IRE; Atwo od et al.[[2009| ), which have the lowest 
background contamination. For each source, we selected only 
photons within a 15° region of interest (Rol) centered around 



' http : //f ermi .gsfc.nasa.gov/ssc/data/analysis/ 


documentation/Cicerone 
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the source itself. To avoid background contamination from the 
bright Earth limb, time intervals where the Earth entered the 
LAT Field of View (FoV) were excluded from the data sample. 
In addition, we excluded observations in which the source un- 
der study was viewed at zenith angles larger than 105°, where 
Earth's atmospheric y-rays increase the background contamina- 
tion. The data were analyzed with a binned maximum likelihood 
technique (Mattox et al. 199 6 ) using the analysis software gtlike 
developed by the LAT tearrj^ A model accounting for the diffuse 
emission and nearby y-ray sources was included in the fit. 

The diffuse foreground, including Galactic interstellar emis- 
sion, extragalactic y-ray emission, and residual CR background, 
was modeled with gZL/em.vO^for the Galactic diffuse emission 
and isotropic Jem_v02 for the extragalactic isotropic emission. 
Each source under study was modeled with a power-law func- 
tion 



Table 13. Summary of y-ray detections with significance TS > 
25. 



dN _ N{Y+\)E^ 
dE ~ _ ^r+i 

nidx mm 



(1) 



where both the normalization factor and the spectral index F 
were allowed to vary in the model fit. The model also includes 
all the sources within a 20° Rol included in Fermi-LAT one-year 
catalog Abdo et al. ( 2010c| l and modeled using power-law func- 
tions. If a source included in the model is a pulsar belonging to 
the Ferm/-LAT pulsar catalog Abdo et al.| ( |2010d[ ), we modeled 
the spectrum with a power-law with exponential cut-off using 
the spectral parameters in the pulsar catalog. 

For the evaluation of the y-ray SEDs, the whole energy range 
from 100 Me V to 300 GeV was divided into two equal logarith- 
mically spaced bins per decade. In each energy bin, the standard 
gtlike binned analysis was applied assuming power-law spectra 
with photon index - -2.0 for all the point sources in the model. 
Assuming that in each energy bin the spectral shape can be ap- 
proximated by a power-law, the flux of the source in all selected 
energy bins was evaluated, requiring in each energy bin a TS 
greater than ten. If the TS is lower than ten, an upper limit (UL) 
was evaluated in that energy bin. Only the statistical errors in 
the fit parameters are shown in the plots. Systematic errors due 
mainly to uncertainties in the LAT effective area derived from 
the on-orbit estimations, are < 5% near 1 GeV, 10% below 
0.1 GeV, and 20% above 10 GeV. 

For each source, we considered the three different integration 
periods for the y-ray data: 

- Simultaneous observations: data accumulated during the pe- 
riod of Planck observation of the source. As the Planck in- 
struments point in slightly different directions and the field of 
view depends on the frequency of observation, a typical ob- 
servation covering all Planck channels takes about one week, 
the exact integration time depending on the position of the 
source. 

- Quasi-simultaneous observations: data integrated over a pe- 
riod of two months centered on the Planck observing period 
of the source. 

- Twenty-seven month Fermi-LAT integration: data integrated 
over a period of 27 months from August 4, 2008 to 
November 4, 2010, i.e., the entire Fermi-hPC[ data set avail- 
able for this paper. 



http : //f ermi . gsf c . nasa . gov/ ssc/data/analy sis/ 

|documentation/Cicerone/Cicerone_Likelihood 

^ http : //f ermi . gsf c . na sa . gov/ ssc/data/access/lat/ 
[BackgroundModels -htmli 





No. of detected 


sources 


Sources 


Sample 


simult. 


2 months 


27 months 


in sample 


Fermi-LAT 


18(36%) 


40 (80%) 


50(100%) 


50 


5w(/f-BAT 


9 (26%) 


12 (35%) 


27 (79%) 


34 


ROSATIKASS 


10 (23%) 


15 (35%) 


24 (56%) 


43 


Radio 


22 (21%) 


38 (37%) 


78 (75%) 


104 



Table 14. Statistics of y-ray detections (TS > 25) in the 27 
month Fermi-LAT data set. 





No. of detected sources 


Sample 


FSRQs 


BL Lacs 


Uncertain 


Fermi-LPC^ 


28 (100%) 


14 (100%) 


8 (100%) 


5wi//^-BAT 


17 (63%) 


7 (100%) 


3 (50%) 


ROSATIKA&S 


8 (53%) 


14 (88%) 


2 (17%) 


Radio 


48 (72%) 


16 (100%) 


9 (64%) 



Tables 13 and 14 give a summary of the y-ray detections 
(TS > 25) in all our samples. The fraction of sources detected by 
Fermi-LPX during the simultaneous integrations is not very large 
and varies from ~ 40% in the Fermi-LAT sample to just ~ 20% 
in the soft X-ray selected sample. We note that even considering 
all the Fermi-LAT data available at the time of writing (27 month 
integration), a sizable fraction of the blazars in the radio and both 
soft and hard X-ray selected samples were not detected. 

Detailed results of the Fermi-hPil analysis are given in 



Tables 15 20 where the observed fluxes or upper limits are 
given in six or three energy bands depending on the source 
brightness. 

Two sources (PKS 0548-322 and NGC 7213) appear as sig- 
nificant y-ray detections in our 27 -month data set, although they 
were n ot included in any of the F ermi-LAT catalogues published 
so far ( |Abdo et al.|2009b||2010cl and[3)- These should therefore 
be considered as new y-ray detections. 



4. The importance of simultaneity 

Blazars are, by definition, highly variable sources. It is there- 
fore important to use simultaneous multi-frequency data to build 
SEDs for comparison with theoretical models. In this section, we 
compare our measurements with data taken from the literature in 
order to derive an estimate of the uncertainties introduced by the 
use of non-simultaneous data in different parts of the spectrum. 

Figure [T] plots the Planck flux density at 44 GHz presented 
in this paper versus the WMAP flux density at 41 GHz from the 
WMAP point source catalogs ( [Bemiett et al.|2003[ [Wright et"ar 
|2009a| l. Some scatter is present, but most of the points lie be- 
tween the two solid lines indicating a factor of two variability. 

Figure [2] plots the X-ray fluxes of the sources observed 
by Swift simultaneously with Planck (see Table 10 1 against 
the X-ray fluxes of the same sources from the BZCAT cata- 
log (Massaro et al. 2009 2010[ ). In this case, a large scatter is 
present, with variations of over a factor of ten. 



http : //heasarc . gsf c . nasa . gov/W3Browse/all/ 
|fermilpsc.html 
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Fig. 1. The Planck 44 GHz flux density of the sources in our sam- 
ple is plotted against the 41 GHz flux density from the WMAP 
five-year catalog (81 sources). The three solid lines represent 
equal flux densities (i.e., no variation) and a factor of two vari- 
ability above or below the equal flux level. Almost all the points 
lie between the factor of two variability lines. 



Figure [3] shows the Fermi-LAT y-ray fluxes of our sources 
measured simultaneously with Planck plotted against their y-ray 
fluxes in the Fermi-LAI IFGL catalog (Abdo et al. 2010c). As 
in the X-ray sample, a scatter with variations larger than a factor 
of ten is observed. 

We conclude that SEDs built with non-simultaneous data 
suff'er from uncertainties in the microwave region that are rel- 
atively modest and generally limited to about a factor of two, 
while the high energy part of the spectrum (X-ray and y-ray) 
is much more affected, with uncertainties caused by flux vari- 
ations of up to a factor of ten or more. The same uncertain- 
ties, of course, apply when searching for correlations in non- 
simultaneous multi-frequency data. 

5. Spectral energy distributions 

We constructed the SEDs of all the blazars in our samples from 
the simultaneous multi-frequency data described above using the 
ASDC SED Builder, an on-line service developed at the ASI 
Science Data Center (ASDCFll jStratta et al.||201l| i. This is a 
WEB-based tool that allows users to build multi-frequency SEDs 
combining data from local catalogs and external services (e.g., 
NED, SDSS, USNO) with the user's own data. The tool con- 
verts observed fluxes or magnitudes into de-reddened fluxes at a 
given frequency using standard recipes that take into account the 
instrument response and assumed average spectral slopes. The 
SED builder can display SEDs both in flux and in luminosity (if 
redshift information is given); it also provides useful tasks such 
as the overlay of templates for blazar host galaxies and nuclear 
optical emission (blue-bump), and allows users to compare the 
SED with models including one or more SSC components. 

The SEDs of all the sources in our samples are shown in 
Figs. 24 41 In these figures, red points represent strictly si- 
multaneous multi-frequency data, green points represent y-ray 



Fig. 2. The Swift X-ray (0. 1-2.4 keV) flux of the sources in 
our sample measured simultaneously with Planck is plotted 
against the 0. 1-2.4 keV flux reported in the BZCAT catalog (83 
sources). The three solid lines represent equal fluxes (i.e., no vari- 
ation) and a factor of two variability above or below the equal 
flux level. Note that several points are outside the factor of two 
variability lines, revealing variability of up to about a factor ten. 




http : //tools . asdc . asi . it/SED/ 



io-« 10-' io-« 10-= 

Fermi flux simultaneous (E > 100 Mev, pli cm-^ s-') 



Fig. 3. The Fermi-LAJ^ y-ray flux of the sources in our samples 
detected during the simultaneous integration with Planck is plot- 
ted against the flux reported in the Ferm/-LAT 1-year catalog. 
The three solid lines represent equal fluxes (i.e., no variation) 
and a factor of two variability above or below the equal flux 
level. Note that several points are outside the factor of two vari- 
ability lines, revealing variability of up to about a factor ten. 



data integrated over a period of two months centered on the 
times of the Swift/ Planck observations, ground-based data taken 
quasi-simultaneously, and P/fl«cA:-ERCSC flux densities, and 
blue points represent y-ray data integrated over the full period 
of 27 months. In the few cases where no Swift simultaneous ob- 
servations could be obtained, we plot only Planck, Fermi-LALT, 
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and XMM-Newton ( [Vagnetti et al.|20l"0 1. To assess the presence 
of a possible thermal component in the X-ray emission, we used 
the relationship given by Grupe et al. ( 2010| l 



Q'UV-X(Radio-quietQSO) - 0- 1 14 log(L2500 a) ~ 1-177, 



(2) 



14.6 14.8 15 

Log (frequency), Hz 



Fig. 4. Top panel: The SDSS template of |Vanden Berk et al.|2001 
for the broad-line and thermal emission from a QSO. Bottom 
panel : The giant elliptical galaxy template of jMannucci et al!] 
120011 See text for details. 



and ground-based data. Two-cr upper limits are indicated by ar- 
rows. 

5.1. Distinguishing tlie non-tliermal/jet-related radiation from 
QSO accretion and liost galaxy emission 



We used the simultaneous SEDs of Figs 24 41] to determine 
some parameters that can constrain the physical mechanisms 
powering blazars. However, before doing so we had to identify 
and separate the radiation that is unrelated to the non-thermal, 
relativistically amplified emission from the jet of the blazars; that 
is, radiation from accretion onto the central black hole and from 
the host galaxy (see, e.g., Perlman et al. 2008 ). To do so, we es- 
timated the contamination by the host galaxy assuming that all 
blazars are hosted by giant elliptical galaxies (e.g., Kotilainen 
let al |1998t|Nilsson et al.|2003j|Le6n-Tavares et al. 2011) with 
absolute magnitude of Mr - -23.7. As [Scarpa et al. (2000,) 



and Urry et al. (2000) demonstrated, this value is within one 
magnitude of the observed values in a sample of over 100 BL 
Lacs observed with the Hubble Space Telescope. For the spec- 



tral shape, we used the elliptical galaxy template of Mannucci 
et al. 2001 (bottom panel of Fig. [4]i, which is based on low 



spectral resolution observations of a number of nearby galax- 
ies in the wavelength range 0.12 -2.4 fim and is a good match to 
the predictions of spectrophotometric models for giant ellipticals 
IMannucci et arilODTj 

The radiation produced by accretion was estimated from the 
composite optical spectrum built by |Vanden Berk et aL] ( |2001| l 
using over 2200 optical spectra of radio-quiet QSOs taken from 
the SDSS database ( [York et al.|[2000| (top panel of Fig |4|, and 
the expected soft X-ray emission of radio quiet AGN from |Grupe] 
let al.120101 

The ratio of optical to soft X-ray light has been known to be a 
function of optical luminosity since the early obervations of the 
Einstein observatory ( Avni & Tananbaum 1986) . More recently, 
this dependence has been confirmed using simultaneously ac- 
quired optical and soft X-ray data from Swift ( Grupe et al. 2010| l 



where ^2500 A 1^ rest-frame luminosity of the thermal emis- 
sion at 2500 A in units of W Hz"' and auv-xcradio-quietoso) is the 
usual slope between the UV (2500 A) and the soft X-ray (2 keV) 
flux in radio quiet QSOs (e.g., [Vagnetti et al.|2010 1. 

Examples of the emission from these components unrelated 
to the jet are shown in Fig.|5j which shows the SEDs of Mkn501 
and Mkn421 where the optical light is dominated by the host 
galaxy. 

Figs. |6 and|7 show the SEDs of BL Lacertae where part of 
the UV light is thought to originate in disk emission (Raiten] 
eraL||2009j), the nearby FSRQ 3C273 (z = 0.158), and the 
high-redshift FSRQs 4C 38.41 (z = 1.814) and IJy 0537-286 
(z = 3.104) where the optical/UV light is heavily, or completely, 
contaminated by radiation coming from accretion onto the cen- 
tral black hole. 

We compared the amplitude of the optical thermal emission 
to non-thermal radiation using the parameter CR-ocThermai), de- 
fined as the spectral slope between the 5 GHz radio flux den- 
sity and the 5000 A optical flux density that can be attributed to 
the blue-bump/disk/thermal emission. This quantity depends on 
both the relativistic amplification factor and the intrinsic ratio of 
non-thermal/jet radiation to disk emission. In the case of FSRQs, 
the optical spectrum displays emission lines by definition, there- 
fore we were able to constrain Q'R-o(Thei-niai) by adjusting the op- 
tical thermal emission to the same level as the data. When the 
thermal blue bump was seen directly in the optical/UV part of 
the SED (e.g., Figs.|6]and|7]), we fit the template of Vanden Berk 



l et al.| ([200T) to the observed data; when the disk emission was 
not obviously visible we set the intensity of the blue-bump tem- 
plate to just below the observed optical/UV emission where the 
broad emission lines had been detected. 

The BL Lacs do not show emission lines in their optical spec- 
trum, so for this class of objects we could only set upper limits on 
ffR-ocThermai)- We did that by assuming that the template for the 
optical thermal emission is at least one order of magnitude below 
the observed data, that is sufficiently low to hamper any broad 
line detection. The estimation of Q'R-o(Thermai) relies strongly on 
the quality of the optical data available, in particular on the si- 
multaneous UVOT data. Therefore we define an "optical data 
quality" flag as follows: 

- 0: no simultaneous data available; 

- 1: poor quality (e.g., only one UVOT filter available); 

- 2: good quality (e.g., two or three UVOT filters available); 

- 3: excellent quality (all UVOT filters available). 

In Figure [8] (upper panel), we show the distribution of 
ffR-O(Theimai) for all the FSRQs with excellent optical data for 
the whole sample and each of our samples independently. The 
distribution of the ffR-ocThermai) upper limits for BL Lac ob- 
jects is shown in the bottom panel. The results suggest a possi- 
ble difference in the Q'R-o(Theimai) distribution between different 
samples. As reported in Table l2T[ when only sources with ex- 
cellent optical data are considered we obtain (aR_o(Thermai)) = 
-0.64 ± 0.05 for the ROSATfRASS sample, to be compared with 
(ffR-ocThei-mai)) = -0.73 ± 0.02 for the Fenni-LALT sample. 

To assess the possible presence of a thermal component in 
the X-ray emission, for each source we compare the predicted 
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10 15 20 25 10 15 20 25 

Log Frequency v (Hz) Log Frequency v (Hz) 



Fig.5. The SEDs of MknSOl (left) and Mkn421 (right) showing the expected emission from the host galaxy (giant elliptical) as an 
orange line. The green lines are the best-fit to the simultaneous non-thermal data using a third degree polynomial function. See text 
for details. 




Log Frequency v (Hz) Log Frequency v (Hz) 



Fig. 6. Left: the SED of BL Lacertae showing the expected emission from the host galaxy (just below the observed non-thermal 
radiation, giant elliptical, orange line) and the blue bump emission (blue line, see also Raiteri et al. 2009| l. Right: the SED of 3C 273 



showing the thermal emission from the blue bump and the expected X-ray emission from accretion including 1, 2, and 3cr bands 
(purple and blue lines) derived from Eq. [2] The vertical parallel lines represent the optical window (4000-10000 A). The green 
lines are the best-fit to the simultaneous non-thermal data using a third degree polynomial function. 

Table 21. Contamination of the X-ray emission from a thermal component and (ffR ocThermai)) estimated for FSRQs in the whole 
sample and for each sample independently, for various thresholds on the quality of the optical data. We have considered as contam- 
inated sources all those with X-ray thermal contamination flag < 2 (see the text for a description of the flags). 









ALL 




Fermi-LAT 




Swift-BKI 




ROSATIKk&S 


Opt Data 


cont/ 


% 


(CR-OCThcrmal)) 


cont/ 


% 


(ttR-OCThermal)) 


cont/ 


% 


(CR-OfThermal)) 


cont/ 


% 


(CR-OCThermal)) 


quality 


total 






total 






total 






total 






ALL 


10/47 


21 


-0.69±0.02 


3/25 


12 


-0.72+0.02 


2/21 


10 


-0.70±0.03 


6/13 


46 


-0.65±0.04 


>2 


10/44 


23 


-0.69±0.02 


3/22 


14 


-0.73±0.02 


2/20 


10 


-0.70±0.03 


6/13 


46 


-0.65±0.04 


3 


9/40 


23 


-0.69±0.02 


3/22 


14 


-0.73±0.02 


1/16 


6 


-0.69±0.04 


6/12 


50 


-0.64±0.05 



thermal emission from accretion with the actual X-ray spectrum, 
and we also include the uncertainties in the parameters of Eq. |2] 
in determining the uncertainties given by 1,2, and 3cr bands. 
We therefore define the following X-ray thermal contamination 
flags: 



- 0: no X-ray data available; 

- 1: X-ray emission mostly or entirely due to accre- 
tion/reflection (data agree with the expectations for accretion 
emission within Icr); 



13 



p. Giommi et al.: Simultaneous Planck, Swift, and Fermi observations of blazars 




Log Frequency v (Hz) Log Frequency v (Hz) 



Fig. 7. The SEDs of 4C 38.41 (z = 1.814; left) and IJy 0537-286 (z = 3.104; right). Simultaneous data ai-e shown in red; non- 
simultaneous literature or archival data are shown in light gray. Note that the UVOT data of these medium and high redshift objects 
matches quite well the QSO template of |Vanden Berk et al.| ( |200 r) that we use to estimate the thermal emission from the blue 
bump (blue line). The emission from the host galaxy (orange line) is very low compared to other components. The observed X-ray 
emission from these sources is more than 3cr above the expected emission from accretion derived from Eq. |2] In both cases, the 
optical light is dominated by radiation from the accretion while the X-rays originate from the non-thermal component. 



Radio sample 
Fermi sampie 
BAT sampie 
RASS sampie 




c 
o 



o 

CO 



T3 
CD 
N 

"c6 

E 



<X3 

d 



CO 

d 



d 



d 



' 1 

Radio sampie 


1 


1 




Fermi sampie 








BAT sampie 








RASS sampie 


























1 






) 2 




4 





'^RO(Thermal) 



X-ray coiitamiiiatioii flag 



Fig. 8. Distribution of a'R-o(Thermai) for all the sources and for 
each blazar sample. Only SEDs with good or excellent optical 
data (flag = 2 or 3 , see text) were used. Black solid histograms: 
radio sample; red dashed histograms: Fermi sample; green dot- 
dashed histograms: Swift BAT sample; blue dotted histograms: 
RASS sample. Distributions for BL Lac objects refer to upper 
limits only. 



Fig. 9. Distribution of the thermal emission flag for all the 
FSRQs with good or excellent optical data flag (see text for 
details). Black solid histograms: radio sample; red dashed 
histograms: Fermi-hPC^ sample; green dot-dashed histograms: 
Swift BAT sample; and blue dotted histograms: RASS sample. 
As the X-ray emission of all BL Lacs was not contaminates, the 
BL Lac data are not plotted here. 



- 2: X-ray emission probably contaminated by the accretion 
component (data agree with with the expectations for accre- 
tion emission within 2cr); 

- 3: X-ray emission mostly of non-thermal origin (data are at 
2-3cr from the expectations for radio quiet QSOs); 

- 4: X-ray emission certainly of non-thermal origin (data are 
more than 3cr away from the expectations for radio quiet 
QSOs). 



Results for all the FSRQs with good or excellent optical 
data are shown in Figure |9j and summarized in Table [2T| for the 
whole sample and each sample independently. We considered as 
contaminated all the sources with X-ray thermal contamination 
flag < 2. There is a large difference between the Fermi-LPsI and 
^wi/f-BAT samples, where < 15% of the sources have a ther- 
mal component in their X-ray emission, and the ROSATIKASS 
sample, where ~ 50% of the sources are contaminated. This 
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demonstrates that the thermal component in the X-ray emission 
of blazars cannot be neglected, even in bright sources. 

5.2. SED parameter estimation 

We used the SEDs of all the objects in our samples to esti- 
mate the values of important physical parameters such as v^^^^, 
yS pryS \ yic ^jj^ ^ic f^yic \ /ggg Table 22 1 taking into 

peak ^ peak^' peak' peak ^ peak^ ^ ^ 

account only the non-thermal radiation and fitting third-degree 
polynomials as described in Abdo et al. ( 2010a[ l (see Figs.|5]|6] 
and|7]for examples). 

For the sources that were not detected by FermZ-LAT even 
in the 27-month integration, we estimated limits of v'^_,,and 



„ic 



peak 



j^F(Vp^^) by constraining the polynomial in the high-energy 
part of the SED with the 27-month ferm/-LAT upper limits, as 
shown in Fig.[TO| 

In some HSP BL Lacs with particularly high v^^^j^ values 
(see, e.g.. Figs. [30] [31] [36] and [38) the Fermi-LAJi data alone 
are insufficient to ensure a good measure of v^^^^, as the spectra 
are still rising at the highest Fermi-LP^ energies and no simulta- 
neous TeV data are available. The v"^ , values for these sources 

peak 

should therefore be considered as lower limits. 




15 20 
Log Frequency v (Hz) 

Fig. 10. The SED of the blazar PKS 0003-066, illustrating the 

estimation of the upper limits to v"-- , and to v^*- , Fiv^'^ , ) by 

peak peak ^ peak^ 

combining the X-ray data with the 27-month Fermi-hKY upper 
limits. 



6. The spectral slope of blazars in the 
radio-microwave region 

While WMAP results are consistent with a single flat spectral in- 
dex in the relatively narrow frequency range 23-94 GHz ( WrighF 
|et al.|2009b |Gold et aL 201 Ij ), the blazar spectrum must steepen 
at frequencies closer to the synchrotron peak. Adding Planck 
data and simultaneous ground-based observations at centimetre 
wavelengths to the WMAP data improves the spectral coverage, 
and allows us to probe the spectral shape of blazars over the 
much wider frequency range ~ 1 GHz to ~ 1 THz. 

We studied the low frequency (LF) and the high frequency 
(HE) regions of the centimetre to sub-millimetre blazar spectra 
separately to search for differences in the spectral index a and 
determine the frequency at which the spectral index changes. We 
fitted the two frequency regions independently with power-laws 



to estimate the spectral indices at both low frequencies (q-lf, for 
y ^ i^Break) and high frequencies (ohp, for v > VBi-eak), assum- 
ing a range of break frequency values, VBieak, from 30 GHz to 
100 GHz. In Fig.[TT] we show the distributions of q;lf and qthf 
for VBieak - 70 GHz and 100 GHz; the blazar spectra steepen 
from (Tlf ~ to ckhf 0.65. 

To verify the robustness of the results we repeated the anal- 
ysis by imposing the two minimum numbers of independent fre- 
quencies needed to perform the fit, namely 3 and 5, and the re- 
sults are consistent (Table [23). 

We also analyzed the different classes of blazars -flat spec- 
trum radio quasars (FSRQ), BL Lac objects (BL Lac), and radio- 
loud AGN of unknown classification (Uncertain Type) - sepa- 
rately, and we have found a difference at a level of about 3cr in 
QfHF, which is ~ -0.7 for FSRQ but ~ -0.5 for BL Lacs. 

The spectral steepening and the difference in ouf are visi- 
ble in Figure 12 where we show ulf versus (vs.) q-hf for the 



sources that meet the requirements for the minimum amount of 
independent data at both low and high frequencies. 

Our results are in complete agreement with the findings of 
IPlanck Collabor ation||201 le| for the radio selected sample, i.e., 
flat spectral index at low frequency and aup 0.6 at high fre- 
quency with a break frequency ~ 70 GHz, suggesting that this 
is a general feature of all blazars regardless of the selection cri- 
teria. There is also general agreement with the results presented 
in IPlanck Co llaboration 20 lid for all the sources in the Planck 
ERCSC catalogs, and the slight discrepancy in the spectral index 
estimated at low frequency can be explained by the fact that, un- 
like Planck Collaboration|201 Id we have included the ground- 
based observations at 5-30 GHz. 
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Fig. 12. (Thf vs. q-lf diagram for vgreak - 70 GHz for different 
blazar classes: flat spectrum radio quasars (FSRQ), BL Lac ob- 
jects (BL Lac), and blazars of uncertain classification. 



7. Searching for correlations between fluxes in 
different energy bands 

We used our simultaneous and average multi-frequency data to 
identify possible correlations between fluxes measured in differ- 
ent energy bands. As some of the blazars in our samples were 
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Table 22. Rest-frame (v^^^) and (v'^^,^) values for different samples and classes. 







Synchrotron peak 


inverse 


Compton peak 


Sample 


Class 


No. of sources 


\^^&\ 1 eakf 1 


No. of sources (yo%(ypeak)) 




ALL 


45 


13.34+0.12 


Al 


22.19+0.10 


Ferrai-LAT 


FSRQs 


24 


13.07+0.07 


23 


11.1\±Q.\\ 




BL Lacs 


13 


14.11+0.32 


11 


22.35+0.24 




ALL 


26 


13.65+0.26 


22 


21.66+0.28 


Swift-BKi 


FSRQs 


18 


13.16+0.09 


16 


21.45+0.17 




BL Lacs 


5 


15.75+0.81 


3 


23.41 + 1.52 




ALL 


25 


14.27+0.31 


20 


22.15+0.38 


ROSAT/RASS 


FSRQs 


11 


13.02+0.17 


11 


21.43+0.30 




BL Lacs 


11 


15.79+0.28 


6 


24.24+0.37 




ALL 


94 


13.20+0.06 


77 


21.99+0.10 


Radio 


FSRQs 


64 


13.08+0.05 


49 


21.99+0.12 




BL Lacs 


14 


13.89+0.27 


13 


22.26+0.20 
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Fig. 11. Radio/microwave low frequency (v < 70 GHz, LF; left side) and high frequency (v > 70 GHz, HF; right side) spectral index 
distributions of FSRQs and BL Lacs in our samples for the case VBreak - 70 GHz. While the distributions of low frequency spectral 

indices are very similar for both types of blazars with (qtlf) 0.1, the distributions of the high frequency slopes indicate that the 

spectra of FSRQs ((qthf) = -0.73 ± 0.04) might be somewhat steeper than those of BL Lacs ((aHp) = -0.5 1 ± 0.07). A Kolmogorov 
Smirnov (KS) test gives a probability of less than 3% that the HF distributions for FSRQs and BL Lacs come from the same parent 
distribution. 



not detected by either Ferm/-LAT or Planck, we estimated the 
significance of the correlations using the ASURV code Revl.2 
( [Lavalley et al. 1992|l, w hich takes account of upper limits as 
described in 'Isobe et al. ( 1986). We search for possible corre- 
lations using fluxes or flux-densities rather than in luminosity- 
luminosity space because this allows us to use all objects in the 
samples, including those with no redshift (and consequently lu- 
minosity) information. 



7 A. Microwave vs. X-ray 



Figure 13 shows the Planck 143 GHz flux density versus the si- 
multaneous Swift-XRT X-ray flux for all the sources where the 
X-rays are expected to be due to the inverse Compton compo- 
nent (that is all LSP blazars) and not significantly contaminated 
by X-ray emission that is unrelated to the jet, such as that pro- 
duced by the accretion process (see Sect. 5.1 for details). Sources 



that were not detected by Planck are plotted as upper limits; all 
the sources were detected in the X-ray band. 



A correlation, although with some scatter, is clearly present. 
The Spearman rank coefficients (p) and the corresponding prob- 
abilities that the observed correlation is the result of chance are 



given in Table 24 Results obtained using Planck flux densities 
at other frequencies from 30 GHz to 217 GHz are similar and not 
shown here. 



7.2. Microwave vs. y-ray 

The relationship between radio or microwave and y-ray fluxes is 
a topic that has been addressed several times in the literature. A 
positive correlation between the radio and y-ray fluxes is gener- 
ally found, though with a large scatter (e.g.,|Kovalev et al" 



2009 



Leon-Tavares et al.|20101 iGiroletti et al.||2()101 Ghirland a et aL| 
2010; M ahony et ai|2010(|Peel et al.|2 010';'Linford et al .|2011t 
Ackermann et al.|2011[|Le6n-Tavares et al. 201 1). However, in 
most cases radio and non-simultaneous y-ray data for sources 
detected in both energy bands are compared. For the first time. 
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Table 23. Radio LF and HF spectral index distributions for different values of the break frequency vgreak and of the minimum 
number of frequency bands considered for the fit. 



Low Frequency High Frequency 



^Break 


n. of frequencies 


Class 


n. of sources 


(q-lf) 


0"LF 


n. of sources 




0"HF 


30 


3 


ALL 


47 


0.00± 0.04 


0.27 


69 


-0.51+0.04 


0.30 


30 


5 


ALL 


46 


0.00± 0.04 


0.27 


57 


-0.49+ 0.03 


0.26 


44 


3 


ALL 


47 


-0.02± 0.04 


0.25 


66 


-0.64+ 0.03 


0.26 


44 


5 


ALL 


47 


-0.02± 0.04 


0.25 


47 


-0.62+ 0.03 


0.22 


70 


3 


ALL 


63 


-0.08± 0.03 


0.28 


66 


-0.67+ 0.03 


0.28 


70 


5 


ALL 


47 


-0.02+ 0.04 


0.25 


35 


-0.65+ 0.04 


0.22 


100 


3 


ALL 


65 


-0.11+0.03 


0.27 


48 


-0.56+ 0.04 


0.26 


100 


5 


ALL 


47 


-0.02± 0.04 


0.24 


15 


-0.52+ 0.05 


0.19 


70 


3 


FSRQ 


39 


-0.1 1± 0.04 


0.28 


40 


-0.73+ 0.04 


0.27 


70 


3 


BLLac 


12 


-0.08± 0.08 


0.29 


14 


-0.51+ 0.07 


0.27 


70 


3 


Unc. Type 


12 


-0.02± 0.08 


0.29 


12 


-0.64+ 0.07 


0.25 



Table 24. Spearman correlation parameters for 143 GHz flux density vs. X-ray flux. 



Class Fenni-hfiH sample Swift-BKT sample ROSAT/RASS sample Radio sample 





No. of 
sources 
(a)/(b) 


P 


P{%) 


No. of 
sources 
(a)/(b) 


P 


P(%) 


No. of 
sources 
(a)/(b) 


P 


P{%) 


No. of 
sources 
(a)/(b) 


P 


P(%) 


All 


35/3 


0.72 


<0.01 


18/4 


0.50 


2.34 


9/10 


0.70 


0.30 


49/3 


0.51 


0.03 


FSRQs 


20/2 


0.75 


0.06 


16/3 


0.58 


1.40 


5/3 


0.89 


1.85 


33/3 


0.57 


0.08 


BL Lacs 


8/1 


0.95 


0.72 


2/0 






1/3 






9/0 


0.28 


42.29 



Number of Planck detections. '*'Number of Planck upper limits. 
* The number of detections is too low to allow a reliable estimation of p. 



E 



► e « ® 



• Radio sample 
O Fermi sample 

★ BAT sample 
□ RASS sample 

FSRQs 

BL Lacs 
Uncertain type 



"100 



1000 

Planck Simultaneous S,, 



10* 



, (mJy) 



Fig. 13. The Planck 143 GHz flux density plotted vs. the simul- 
taneous Swift XRT X-ray flux for all sources where the X-ray 
flux is expected to be due to the inverse Compton component. 
Blazars with significant X-ray contamination due to accretion 
have been excluded from the plot. 



we present simultaneous microwave and y-ray data and take into 
account upper limits. 

The top panel of Figure [14] shows the simultaneous Fermi- 
LAT y-ray flux plotted versus the Planck flux density at 
143 GHz. Sources with TS < 25 in the Ferm/-LAT data and 
below 4cr in Planck maps are shown as upper limits. The 
plot shows a clear trend but from Table 25 which gives the 
Spearman's p correlation parameter and the probability P that 
the observed level of correlation is caused by chance, we see 



that the level of significance of the coiTelation is never very high 
{P is of the order of a few percent and never lower than 0.05%), 
especially in the hard X-ray and y-ray selected samples and for 
BL Lacs. This result is partly due to the large number of upper 
limits to the simultaneous y-ray flux, which was estimated typi- 
cally over a period of one week. To improve the statistics, in the 
middle panel of Fig. [14] we show the same plot using the Fermi- 
LAT y-ray flux integrated over a two-month period centered on 
the time of the Planck observations. From Table [25] we see that, 
although the significance of the correlation in the various sam- 
ples and blazar classes substantially increases, it is never very 
high and the scatter remains. 

This result may be due to the limited flux dynamic range 
in the y-ray band for the FermZ-LAT sample and to the small 
number of objects (in the case of BL Lacs). However, a weak 
coiTelation could occur if the microwave flux density represents 
the superposition of multiple synchrotron components, while the 
simultaneous y-ray flux represents the emission from a single 
dominant component that may be active for only a short time. 
For this reason, in the bottom panel of Fig. 14 we plot the Planck 
ERCSC flux densities (which in most cases represent the flux 
density averaged over more than one Planck observation) versus 
the ferm/-LAT flux averaged over the entire 27-month period. 
In this case the correlation is highly significant, although, the 
many upper limits in both energy bands clearly imply that the 
dispersion is very large. 



7.3. X-ray vs. y-ray 

Figure [15] plots the 2-month y-ray flux versus the simultaneous 
X-ray flux for all sources observed by Swift that do not show sig- 



17 



p. Giommi et al.: Simultaneous Planck, Swift, and Fermi observations of blazars 
Table 25. Microwave (143 GHz) vs. y-ray flux correlation parameters. 



y-ray Class Fermi-LKi sample Swift-BPS sample ROSATfRASS sample Radio sample 



integration 




No. of 


P 


P{%) 


No. of 


P 


P{%) 


No. of 


P 


P{%) 


No. of 


P 


P{%) 


period 




sources 
(a)/(b)/(c) 






sources 
(a)/(b)/(c) 






sources 
(a)/(b)/(c) 






sources 
(a)/(b)/(c) 








All 


16/5/32 


0.22 


124 


8/6/20 


0.42 


2.49 


6/20/25 


0.42 


1.42 


22/2/71 


0.35 


0.07 


Simult. 


FSRQs 


10/3/16 


0.23 


234 


5/2/12 


0.39 


9.92 


4/3/6 


0.65 


3.97 


15/1/46 


0.33 


1.19 




BL Lacs 


5/2/9 


0.55 


4.76 


2/2/2 






2/11/4 






4/0/10 


0.38 


17.6 




All 


36/5/9 


0.29 


4.18 


11/7/20 


0.61 


0.07 


8/21/18 


0.31 


6.00 


37/3/56 


0.40 


0.01 


2 months 


FSRQs 


19/3/6 


0.17 


38.7 


6/3/15 


0.59 


0.89 


4/4/9 


0.51 


6.83 


23/2/42 


0.36 


0.39 




BL Lacs 


11/2/1 


0.72 


0.92 


3/3/2 






3/11/5 






10/0/6 


0.42 


10.2 




All 


47/3/0 


0.47 


0.10 


22/16/7 


0.67 


<0.01 


12/23/18 


0.40 


0.99 


75/4/16 


0.48 


<0.01 


27 months 


FSRQs 


25/3/0 


0.48 


1.27 


16/3/4 


0.54 


1.59 


8/3/7 


0.85 


0.15 


47/2/19 


0.48 


0.01 




BL Lacs 


14/0/0 


0.72 


0.91 


3/4/0 






2/14/2 






16/0/0 


0.55 


3.22 



Number of sources detected both by Planck and Ferra;-LAT. "''Number of Planck upper limits. '^'Number of Fermi-hAT upper limits. 
* The number of detections is too low to allow a reliable estimation of p. 
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Fig. 14. Top panel: The Planck flux density at 143 GHz is plot- 
ted against the simultaneous Fermi-LAI flux. Middle panel : The 
Planck flux density at 143 GHz is plotted against the Fermi-LAI 
flux integrated over the 2-month peiiod centered on the Planck- 
Swift observations. Bottom panel: The Planck ERCSC flux den- 
sity at 143 GHz (flux averaged over more than one Planck sur- 
vey) is plotted versus the y-ray fluxes averaged over the entire 
27 month Fermi-LKI observing period. 



natures of thermal contamination in the X-ray spectrum. Open 
red circles represent HSP sources, where the X-ray flux is due to 
the tail of the synchrotron emission, while black filled circles are 



LSP and ISP sources for which the X-ray flux is related to the 
inverse Compton radiation. This distinction is needed in order to 
properly compare sources where the emission is produced by the 
same mechanism. We therefore compute the y-ray vs. X-ray cor- 
relation coefficient only for LSP and ISP sources. The Spearman 
rank test shows moderate evidence of a correlation that is less 
significant in the longer integrations: P - 2.55% for simultane- 
ous data, P - 6.9% for a 2-month integration, or P = 17.4% for 
a 27 -month integration. 

Figure [16] plots the power-law spectral index in the 
Swift XRT energy band (0.3-10 keV) versus the spectral slope 
in the Fermi-LAI y-ray band ( 1 00 Me V-1 00 Ge V) derived using 
the 2-month y-ray data. A significant correlation is present, con- 
firming the results of Abdo et al.^(.2010a) . The top right corner of 
Fig. 16 shows the same plot built with simultaneous y-ray data; 



although the number of objects is smaller, the coiTelation is still 
present. 



8. Discussion and conclusions 

We have collected simultaneous Planck, Swift, Fermi-LAI, and 
ground-based multi-frequency data for 105 blazars included in 
three statistically well-defined samples characterized by flux 
limits in the soft X-ray (0.1-2.4 keV, ROSAT), hard X-ray (15- 
150keV, Swift-B AI), and y-ray (E > 100 Me V, Fermi-LAI) en- 
ergy bands, with the addition of a cut to the radio 5 GHz flux 
density to ensure a high probability of detection by Planck. This 
study complements a similar study of 104 radio-bright AGN 
(hicHz > IJy) ( |Planck Collaboration]|201 le| l. Altogether, the 
four samples contains a total of 175 distinct objects. The acqui- 
sition of this unprecedented multi-frequency/multi-satellite data 
set was possible thanks to cooperation between the Planck, Swift, 
and Ferm/-LAT teams, who agreed to share data and organize 
an extensive program of multi-frequency observations involving 
over 160 Swift ToO pointings scheduled when the blazars were 
within the field of view of the Planck instruments. 

We have used this unique multi-frequency dataset to build 
well-sampled, simultaneous SEDs of all the blazars included in 
our high-energy selected samples. This collection of SEDs is 
an improvement over previous compilations (e.g., |Abdo et ^ 
2010a) because: a) the SEDs presented here are strictly simul- 
taneous, while in Abdo et al. (2010ai the multi-frequency data 
were collected over a period of up to nine months centered on 
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Fig. 15. The Swift XRT flux is plotted against the two-month 
Ferm/-LAT flux for the sources included in the three high-energy 
flux-limited samples. Open red circles represent HSP sources, 
i.e., high synchrotron-peaked BL Lacs where the X-ray flux is 
dominated by synchrotron radiation. Filled black circles repre- 
sent LSP sources, i.e., blazars with low synchrotron peak, for 
which the X-ray flux is dominated by inverse Compton radia- 
tion. 




as the light from the blazar's host galaxy and the radiation pro- 
duced by the accretion onto the central black hole, which often 
contaminate the non-thermal blazar spectrum in the optical, UV, 
and X-ray bands. 



Our findings are broadly consistent with those of jAbdo et al. 
( 2010a[ l. However, our analyses of larger samples selected in dif- 



-3 -2 
Photon index (y-ray band) 



ferent parts of the electromagnetic spectrum, wider wavelength 
coverages, different level of simultaneity, and the ability to sep- 
arate the emission components, have allowed us to make signif- 
icant progress in several areas. The use of four widely different 
samples has allowed us to investigate the consequences of se- 
lection effects on the estimation of critical parameters such as 
(Vp^^j,), average Compton dominance, especially for the case of 
BL Lacs. 

Some of our sources have been observed simultaneously by 
Swift, Planck, and Fenni-LALT during more than one Planck sur- 
vey. In these cases, we have presented only the data collected 
during the first observation. Multiple simultaneous observations 
of a subset of our blazars and detailed model fitting of the SEDs 
will be the subject of future papers. The main results of this work 
are discussed below. 



8.1. Fermi-LAT detection statistics and the effects of 
variability 

The percentage of Fenni-LALT detected sources during the simul- 
taneous integrations, typically lasting about one week, ranges 
from < 40% in the Ferm;-LAT sample to 20-25% in the radio 
and X-ray selected samples (see Table [T3] l. When 2-month in- 
tegrations centered on the Planck-Swift observations are consid- 
ered, these percentages grow to 80% in the Fenni-LALT sample 
and to ~ 35% in the other samples. However, even when using 
data from the 27-month Ferm/-LAT integrations available at the 
time of writing, many of the blazars belonging to the radio, and 
both soft and hard X-ray selected samples remain undetected. 

We note that the detection rate is quite different for FSRQs 
and BL Lacs: if we exclude the Fenni-LAJi sample where all 
the objects have been detected (by definition), the percentage of 
detections is 95% for BL Lacs and only ~60% for FSRQs, with 
values ranging from 72% in the radio sample to just 53% in the 
RASS sample (see Table [l4|. 

A comparison of our simultaneous data with published and 
archival measurements shows that the use of non-simultaneous 
data in the SED of blazars typically introduces a scatter of about 
a factor of two in the microwave band, and a factor of up to ten 
or more in the X-ray and y-ray bands. 



Fig. 16. The power-law spectral index in the Swift XRT energy 
band (0.3-lOkeV) is plotted versus the slope in the Fenni-LALT 
y-ray band (lOOMeV-lOOGeV) derived using Fenni-LALT data 
integrated over a period of 2 months. The same plot built with si- 
multaneous y-ray data is shown in the inset at the upper right. A 
clear anti-correlation is present; the Spearman test gives a prob- 
ability of less than 0.01% that the correlation is due to chance, 
even in the simultaneous data. 



the first three months of Fermi-LAI operations; b) the sources 
were selected according to different statistical criteria allowing 
us to probe the blazar parameter space from widely different 
viewpoints; and c) we took care to identify and separate radi- 
ation components unrelated to the emission from the jet, such 



8.2. The spectral slope in the radio-sub-millimetre region 

We confirm that the energy spectrum of blazars in the radio- 
microwave spectral region is quite flat, with an average slope 
of {a} ~ (/(y) oc V ") up to about 70 GHz, above which it 
steepens to (a) 0.6. This behaviour is very similar to that ob- 
served in the sample of radio bright blazars considered by |Planck| 
Collaboration|( 20^1 le[) ( see also Tucci et al.|201U confirming the 
findings of Abdo et al. ( 2010a[ ) that the radio to microwave part 
of the spectrum is approximately the same in all blazars (FSRQs 
and BL Lacs) independently of the selection band. However, the 
spectral slope of BL Lacs above ~70 GHz is flatter than that of 
FSRQs with (q-hf) = -0.51 ± 0.07 for BL Lacs compared to 
(ffjjp) = -0.73 ± 0.04 for FSRQs (see Table [23] and Fig.fTT). A 



KS test, performed on the subsamples of FSRQs and BL Lacs 
with radio flux density larger than 1 Jy, gives a probability of 
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less than 3% that the two samples come from the same parent 
population. This difference in the high frequency spectral index 
may reflect that the radio-submm part of the spectrum is closer 
to Vpg^j^ in LBL than in HSP sources. 



8.3. Synchrotron self-absorption 

We searched for signatures of synchrotron self-absorption, 
which in simple homogeneous SSC models is expected to cause 
strong spectral flattening below ~ 100 GHz, but we found no ev- 
idence of any common behaviour; indeed, the average spectrum 
in that region steepens instead of flattening. Possible cases where 
some evidence of synchrotron self-absorption may be present are 
PKS 0454-234, PKS 0521-36 (Fig. 
PKS 1127-145 (Fig. [31]), and 3C45. 



27 1, 8409174-44 (Fig. 29 1, 



.3 (Fig. 40 1 



8.4. Non-thermal versus disk radiation 

In several blazars, the optical/UV light is contaminated signifi- 
cantly by thermal/disk radiation (known as the blue bump, see 
Figs. |6]and Tj, while the soft X-ray flux is contaminated by ra- 
diation produced in the accretion process in approximately 25% 
of the blazars in our samples (see Fig. |9]l. In some of the clos- 
est sources, the optical light is instead either contaminated or 
dominated by the emission from the host galaxy. Ignoring this 
contamination may cause an overestimate of the position of both 

^peak^'^'^ ^peak^y '^^^ '^^ mOrC. 

We investigated the relationship between the radiation pro- 
duced by accretion and the jet, which co-exist in most FSRQs, 
using the parameter aR-o(Theimai), defined as the spectral slope 
between the 5 GHz radio (non-thermal) flux and the 5000 A op- 
tical flux that can be attributed to the blue-bump/disk emis- 
sion. In the blazar paradigm, this quantity depends on both the 
amount of relativistic amplification (which only affects the non- 
thermal radiation from the jet) and the intrinsic ratio of non- 
thermal/jet radiation to disk emission. We estimated the value 
of CR-ocThermai) in all sourccs for which we had optical data of 
good quality. Since BL Lacs do not display broad lines, only 
lower limits to Q;R-o{Thennai) can be derived. 

Figure M shows the distribution of CR-ocThermai) for all our 
samples, which range from ~0.4 to just above 1.0 and have 
peak values between 0.6 and 0.8. The distributions are all sim- 
ilar, with the largest diff'erence being between the Fermi-LAT 
and the RASS samples (a KS test gives a probability of 3.6% 
that the two distributions originate from the same parent popu- 
lation), possibly reflecting differences in either the amplification 
factor or the ratio of accretion to jet emission between y-ray and 
soft X-ray selected blazars. We note that the distribution of the 
Q'R-o(Thermai) limits for BL Lacs is consistent with intrinsic values 
of ffR-ocThermai) for BL Lacs that are within the range of values 
observed in FSRQs. 



8.5. The distribution of rest-frame peak energies 

The distribution of rest-frame synchrotron peak energies (Vp^^j^) 
of FSRQs is very similar in all our samples with a strong peak 
at a: 10'^'^ Hz, an average of (Vpejjj)lO'"' ''''^ ' Hz, and a disper- 
sion of only ~ 0.5 dex. However, for BL Lacs the value of {Vp^^j^) 
is at least one order of magnitude larger than that of FSRQs, 
the exact value depending considerably on the selection method 
(see Table [22] and Fig. 17 1. Since all the sources that are be- 
low the radio flux density cut in the RASS and BAT samples 



are ISP or HSP blazars, their inclusion would increase the dif- 
ference between the , distributions. The distributions of v"-- , 

peak peak 

for FSRQs and BL Lacs also differ, but not as much as those of 



peak ^^^^ ^^r majority of the sources in all the samples 



(both FSRQs and BL Lacs) peak between 10^' Hz and 10^^ Hz, 
with a few extreme HSP BL Lacs reaching ~ 10^^ Hz. 



Radio sample 
Fermi sample 
BAT sample 
RASS sample 




Fig. 17. The rest-frame synchrotron Vpeak distributions of FSRQs 
and BL Lacs in different samples. Black solid histograms: radio 
sample; red dashed histograms: ferm/-LAT sample; green dot- 
dashed histograms: Swift BAT sample; blue dotted histograms: 
RASS sample. 



Radio sample 
Fermi sample 
BAT sample 
RASS sample 




Fig. 18. The rest-frame inverse Compton Vpeak distributions 
of FSRQs and BL Lacs in different samples. Black solid his- 
tograms: radio sample; red dashed histograms: ferm/-LAT sam- 
ple; green dot-dashed histograms: Swift BAT sample; blue dotted 
histograms: RASS sample. 
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8.6. Correlations between fluxes and other blazar 
parameters 

Despite the strict simultaneity of our data, plots of fluxes in dif- 
ferent spectral regions (microwave vs. X-ray, microwave vs. y- 
ray, and X-ray vs. y-ray) still have a large scatter (see Figs. [13] 



tions between the radio and y-rays have been reported (e.g. 
Kovalev et al. 2009; Leon-Tavares et al. 2010; Giroletti et al. 



2010 . Ghirlanda et al.,,2010, ,Mahony et al ,,2010, ,Peel et al. 
20Tol ILinford et al.||201 ![ [Ackermann et al.||2011| l' The differ 



ence might be due to the different synchrotron peak energies of 
the objects in the samples. Even in simple SSC scenarios, this 
introduces a scatter in the correlation between the fluxes (e.g., 
for the same radio flux, an object with higher Vp^^j, is expected 

to produce more y-rays than one with a smaller Vp^^j^). The large 



scatter present in Fig. 14 could also imply that y-ray emission 



is due to components not always directly related to radiation in 
other energy bands, e.g., multiple SSC components (see also 
Abdo et aL||20lOa i. A good correlation is however present be- 
tween the X-ray and y-ray spectral slopes (see Fig.[T6]l. 

We confirm the correlations between the Ferm/-LAT s pectral 
index and the SED peak energies , and v^"-- , found by 

^ ^ peak peak 



Abdo 



et al. (2010a I. As an illustration of the agreement. Fig. 19 plots 
the Fermi-LAT spectral slope estimated using the full 27-month 
data set as a function of v^^^^. The gray points represent the y- 
ray spectral slopes estimated using the quasi-simultaneous two- 
month dataset. These points, plotted without the much larger er- 
ror bars to avoid confusion, clearly cluster around the 27-month 
data, confirming the correlation. The solid line represents the 
best-fit obtained by | Abdo et"aL| ( |2010a| l. 



8.7. Comparison with the expectation of simple SSC models 
As discussed in 



Abdo et al. 



( 2010a[ l, simple SSC models pre- 
dict that in the Thomson regime the peak frequency of the syn- 
chrotron (Vpg^j.) and inverse Compton (Vp^^j^) components are re- 
lated by 



14 and 15 1. This is somewhat surprising, as positive correla- ic 



peak 



peak 



(^peak) 



(3) 



where J^^^ is the Lorentz factor of the electrons radiating at the 
peak energy. This is related to the observed peak frequency of 
the observed photon spectrum by 



^peak ^ 



1/2 



peak 



(4) 



where Vp^^|^ is the synchrotron peak frequency in the rest-frame 
of the emitting region, B is the magnetic field, and 6 is the usual 
Doppler factor (e.g. Urry & Padovani|1995| l. 

In objects where Vp^^j^is higher than ^ lO'^ Hz, the Thomson 
approximation is no longer valid and the inverse Compton 
scattering occurs under the Klein-Nishina (KN) regime. Using 
Monte Carlo simulations, |Abdo et al.| ( |2010a| l estimated the area 
covered by SSC models in the plane log(Vp^^j,)-log(yp^^^^). This 



area is delimited by the solid contour lines shown in Fig. 20 



where we plot the log(ypg^) of our sources, calculated from the 
observed values of Vp^^j^ and v^^^|^ using eq. jsj versus the rest- 
frame log(v^^^j^). 
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Fig. 19. The y-ray spectral index, estimated from the entire 27- 
month Fermi-LAI data set of all sources in our samples is plot- 
ted against log(Vp^^). Gray points, plotted without statistical er- 
rors to avoid confusion, represent the Ferm/-LAT spectral slopes 
estimated using the 2-month integrations. The solid line is the 
best-fit given by [Abdo etaLl(|2010a i 
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Fig. 20. log(y^^^), obtained from Eq. jsj for all the objects in 
our samples for which we could obtain , and v"- , is plot- 

^ peak peak ^ 

ted against log(Vp^^|^) in the rest-frame of the blazars. The two 
black lines delimit the area predicted by simple homogeneous 
SSC models obtained through extensive Monte Carlo simula- 
tions (see |Abdo et al.|2010a| for details). 



As in the case of the bright y-ray blazars considered by Abdo 
|et al.| ( [2010a| l, only a few objects are inside or close to the SSC 
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Fig. 21. The logarithm of the Compton dominance is plotted as 
a function of lo2(v^ , ) for all sources for which , and v"^ , 

peak^ peak peak 

could be reliably determined. 



area, implying that simple SSC models cannot explain the SED 
of many blazars in our samples. This conclusion is supported by 
the lack of a strong correlation between radio and y-ray fluxes. 
However, the Tpf^^ of blazars that were not detected by 

^ic 

peak 



Fermi-LAI and for which we could only infer a limit to 
(from 30% to 40% of the FSRQs in the radio and X-ray selected 
samples; see Fig. 10 for one example) are plotted as upper lim- 
its in Fig. 20 many of these limits are close to or inside the SSC 



area and therefore the SEDs of these objects are likely consistent 
with simple SSC emission. 

8.8. The Compton dominance of blazars 

The Compton dominance (CD, defined as the ratio of the inverse 
Compton to synchrotron peak luminosities) is a crucial parame- 
ter for the study of blazar physics, as it is strictly related to the 
location of the maximum power output in the energy spectrum 
of a blazar. 



Figure 21 plots the CD values, estimated from our SEDs, as a 
function of v^^^^, showing that log(CD) ranges from about -0.5 
to about 2. The larger values are always associated with LSP 
objects, while HSP sources always have values of log(CD) lower 
than ^ 0.5. In this figure, the two blazar subclasses appear to 
be quite different, with BL Lacs having significantly smaller CD 
values, even when their v^^^j, values are equal to those of FSRQs. 
To better understand this difference, in Fig. 22 we plot the CD 
distributions of FSRQs and BL Lacs for different samples and 

, intervals. 

peak 

The FSRQs included in the FermZ-LAT sample, which are 
y-ray bright by definition, show a CD distribution peaking at 
large values. We note that in this sample we also applied a ra- 
dio flux-density limit of 1 Jy, hence the sources below the radio 
cut must be on average more Compton-dominated than those in 
our sample. This implies that the distribution of purely y-ray se- 
lected blazars must be even more strongly peaked at high CD 
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Fig. 22. Distributions of the Compton dominance for FSRQs and 
BL Lacs of the LSP and HSP type. A significant fraction of 
FSRQs in the radio, RASS, and BAT samples have not been de- 
tected in the y-ray band and therefore only limits to the Compton 
dominance (shown as a dashed histogram; only the radio sample 
to avoid confusion) can be calculated. The large difference be- 
tween the CD distribution of FSRQs in the Ferm/-LAT and the 
other samples illustrates the strong bias that y-ray selection in- 
duces. 



values than that of the Fermi-LA^ sample. Considering instead 
FSRQs selected in the radio and the X-ray bands, we get quite a 
different picture, with a broader distribution extending to values 
of less than 1. Moreover, about 30^5% of FSRQs in the radio, 
soft X-ray, and hard X-ray selected samples are not detected by 
Fenni-LALT and therefore they must populate the part of the CD 
distribution with low CD values. This is shown by the dotted red 
histogram, which also includes upper limits to the CD estimated 
as the ratio of the upper limit to v"-- , F(v^'^ , ) and , F(v^ , ) 

peak peak npak ^ np"i-' 



peak ^ peak"^ 



where limits to v^"-- , and v^'^ vF(v^'~ , ) are obtained by fitting the 

peak peak ^ peak^ <^ 

X-ray data together with the 27-month Fenni-LAT upper limits 



as shown in Fig. 10 



8.9. The blazar sequence 



The top panel of Fig. 23 plots the logarithm of the bolometric 
power, represented by the sum of the synchrotron and inverse 

Compton peak luminosities [LboI ~ ^peak^(^peak) + ^peak^(^peak)] 

as a function of log(Vp^^|^) for all sources in the four samples 
considered in this paper and Planck Collaboration (20 lie) for 
which an estimate of , and the bolometric luminosity was 

peak 

possible. We use this plot to test the relationship known as the 
Blazar Sequence, which is the strong anti -correlation between 



Fossatietal. (1998i 
^ject of hvely debate 
2003; Caccianiga & 
2006, Padovani 2007, Ghisenim 

i I II If I 



bolometric luminosity and v^^^i^claimed by 
and |Ghisellini et aL ](|1998[) that remains a su 
(e.g., Giommi et al. 1999; Padovani et al 
Marcha^ 2004j Nieppola et al, 
|& Tavecchio|2b08| ). 

We stress that to robustly test for the existence of such a re- 
lationship it is mandatory to use samples that are unbiased, that 
is selected in such a way that no particular part of the L^d-v^^^y^ 
diagram is more likely than others to be selected. Although our 
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Fig. 23. Top panel: The bolometric luminosity (represented by 
the sum of synchrotron and inverse Compton peak luminosities) 
is plotted against v^^^j^. The non-thermal optical light of blazars 
with no strong emission lines above the dashed line would be 
bright enough to swamp the emission from the host galaxy mak- 
ing the source appear featureless thus hampering any redshift 
measurement. The lower limits, representing BL Lacs with no 
known redshift, are estimated assuming that the non-thermal 
light is ten times brighter than that of a giant elliptical host 
galaxy (see text for details). The bottom panel shows the same 
plot for the subsample of sources satisfying the selection criteria 
of |Fossati et aL] ( |T998] l. 



samples are statistically well-defined, they are not unbiased from 
this viewpoint for the following reasons: 

and 



18 



the distribution of 



peak 



a) As shown in Figs. HTI ; 
strongly depends on the selection rneffiod; hence, for the various 
samples, we obtain different samplings of the parameter v^^^i^. 

b) The area above the dashed Une, which represents the lumi- 
nosity above which the non-thermal emission of a blazar com- 
pletely dominates the observed optical flux{^ cannot be popu- 
lated by blazars with no emission lines, such as BL Lacs, as in 
this case they would appear completely featureless and therefore 
their redshift could not be measured. In this respect, we note 
that over 40% of the BL Lacs in the BZCAT catalog, and an 



We assume a non-thermal luminosity one order of magnitude higher 
than that of the host giant-elliptical galaxy of luminosity equal to that 
found by |Scarpa et al.| ( |2000^ and |Urry et"ari ( |2000^ and verified by us 
to fit our SEDs. 



even larger fraction of the Ferm;-LAT detected BL Lacs, still 
lack any redshift measurement (Mass aro et al.|2010 Shaw et al. 
|2009[[201 0). All the BL Lacs in our samples for which redshifts 
are unknown are plotted in Fig.|23]as lower limits estimated as- 
suming that their non-thermal light is ten times brighter than the 
optical light of the host galaxy. Some BL Lacs of known lu- 
minosity (red points in Fig. [23] l are above the red line because 
their redshift was measured from emission lines with rest-frame 
equivalent widths below the 5A limit; they might be objects with 
properties in-between those of BL Lacs and FSRQs (see. e.g., 
Ghis ellini et al.|201 1| ). 

c) The different radio flux-density cuts applied to our sam- 
ples imply that different subsamples probe different parts of the 
radio luminosity function. The radio and y-ray selected samples, 
which are defined using a high radio flux limit {S >\ Jy), probe 
the high-luminosity end of the radio luminosity function, which, 
for , - 10'^ Hz and the observed redshift and CD distribu- 

peak 

tions, translates into LboI ^ 10^^ erg s"'. The soft and hard X-ray 
selected samples have a radio flux density cut of S - 0.1-0.2 Jy, 
or about one order of magnitude fainter than that of the radio 
and Fermi-LAT samples, hence might include significantly less 
powerful sources, as faint as bolometric luminosities of the order 
of 1 044 ergs-'. 

If we remove the Vp^^^ dependence of the selection method 
by considering only FSRQs (which have the same v^^^j, distri- 
butions as all samples), we see that the luminosity values span 
five orders of magnitude from ~ lO^^ to ~ 10^^ erg s ' and show 
no trend with v^^_^^, which only ranges between ~ lO'^^ and 

~ 10'4 Hz. No obvious correlation is present in each sample sep- 
arately or in the union of the four samples. The L-shaped distri- 
bution that is apparent i n Fig. [23] if lower limits are ignored, is 

similar to that found by Meyer et al. 1(1201 lb who estimated both 

„s 



the Vpg^j^ and peak luminosities of a large sample of blazars us- 
ing non-simultaneous multi-frequency data. These authors, how- 
ever, instead of considering lower limits to the peak luminosities 
of blazars with unknown redshifts, assumed a luminosity corre- 
sponding to the redshift that the host galaxy would have for the 
observed blazar optical magnitude. |Meyer et al. ( 2011| l argued 
that the strong correlation predicted by the blazar sequence turns 
into a blazar envelope when partly misaligned blazars are in- 
cluded in the samples. However, Giommi et al. (2012 ), by means 
of detailed Monte Carlo simulations, showed that this envelope, 
or L-shaped distribution, is expected when blazars with no red- 
shift measurements are not properly taken into account. 

Finally, we consider the subsample of sources that satisfy 
the same conditions as in [Fossati et al.| ( |1998| l, that is 5 5 ghz > 
2Jy for FSRQs, ^sghz > 1 Jy for radio-selected BL Lacs, no 
restrictions for X-ray selected BL Lacs, and the exclusion of all 
BL Lacs with no redshift information. This case is illustrated in 



Fig. 23 (bottom panel), which shows a trend very similar to that 
presented in Fossati et al. ( 1998| l. 



Taking take into account all of the above, we conclude that 
our data do not show a correlation of the type predicted by 
Fossati et al. ( 19981), owing to the presence of low luminosity 



LSP objects and the difficulty in measuring the redshifts of likely 
high-luminosity HSP sources. That such a correlation becomes 
evident when the objects are selected by the criteria of Fossati] 
|et al. ( 1998| l supports the hypothesis that the correlation might 
be the result of a selection effect. 
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8.10. Selection effects and sample composition 

Our decision to select flux-limited samples for widely different 
parts of the electromagnetic spectrum (radio, soft X-ray, hard 
X-ray, and y-ray) has allowed us to demonstrate the strong se- 
lection biases that can affect important physical parameters, such 
as the peak energy of both the synchrotron and inverse Compton 
components (see Figs. [17] and 18 1 and the Compton dominance 
(see Fig. 22 1. Since FSRQs and BL Lacs have significantly dif- 



ferent Vpg^j. distributions, these selection biases also strongly 
affect the composition of the samples in terms of the relative 
abundances of blazar subclasses (FSRQs vs. BL Lacs, LSPs vs. 
HSPs), as is apparent from Table |4] 

Radio-selected samples include sources that are bright 
in the radio band. If there is no correlation between radio 
ffux/luminosity and other parameters such as v^^^i^ and Compton 
dominance, this is the best selection for measuring the distri- 
butions of these important physical parameters. If instead there 
is a strong correlation between radio luminosity and v^^^y_, then 

the distribution of v^^^^^ should strongly depend on the radio flux 
limit. 

X-ray selection favors high Vp^^i^ (and consequently high 
^eak-* sources, which are much brighter at X-ray frequencies 
than low v^^^j^ for the same radio flux. X-ray flux-limited samples 
are therefore much richer in high v^^^i^ BL Lacs (HBLs or HSP 
sources) than radio-selected samples. This selection effect has 
been known since the first soft X-ray surveys became available. 

Selection in the y-ray band favors bright y-ray objects and 
therefore highly Compton-dominated sources. ferm/-LAT TS- 
limited samples contain more sources with flat y-ray spectral 
slopes or high Vp^^j,sources. This explains the overabundance of 
HSP blazars (only BL Lacs) and high CD blazars (only FSRQs) 
in Fermi-LAiT catalogs. 
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Fig. 24. The SED of IIIZW2 (JOOlO+1058, top left), S5 0014+813 (J0017+8135, top right), lES 0033+595 (J0035+5950, middle 
left), Mkn348 (J0048+3157, middle right), IJy 01 18-272 (J0120-2701, bottom left), and S40133+47 (JOl 36+4751, bottom right). 
Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck ERCSC and non- 
simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in blue; literature or 
archival data are shown in light gray. 
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Fig.25. The SED of PKS0202-17 (J0204-1701, top left), PKS0208-512 (J0210-5101, top right), GB6J 0214+5145 
(J0214+5144, middle left), PKS 0215+015 (J0217+0144, middle right), Uy 0212+735 (J0217+7349, bottom left), and 
IJy 0218+357 (J0221+3556, bottom right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data inte- 
grated over 2 months, Planck ERCSC and non-simultaneous ground based observations are shown in green; Fermi data integrated 
over 27 months are shown in blue; literature or archival data are shown in light gray. 
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Fig. 26. The SED of 4C 28.07 (J0237+2848, top left), PKS 0235+164 (J0238+1636, top right), NGC 1275 (J0319+4130, middle 
left), PKS 0332-403 (J0334-4008, middle right), NRAO 140 (J0336+3218, bottom left) and PKS 0420-01 (J0423-0120, bottom 
right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck ERCSC and 
non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in blue; literature 
or archival data are shown in light gray. 
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Fig. 27. The SED of PKS 0426-380 (J0428-3756, top left), 3C 120 (J0433+0521, top right), PKS 0454-234 (J0457-2324, middle 
left), PKS 0521-36 (J0522-3627, middle right), PKS 0528+134 (J0530+1331, bottom left), and PKS 0537-441 (J0538-4405, 
bottom right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck 
ERCSC and non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in 
blue; literature or archival data are shown in light gray. 



29 



p. Giommi et al.: Simultaneous Planck, Swift, and Fermi observations of blazars 
J0539-2839 - 1Jy 0537-286 J0550-3216 - PKS0548-322 




1 1 1 1 1 

i 
i 






• 




i 






444- 




1 1 









15 20 
Log Frequency v (Hz) 

J0623-6436 - I R AS- L 06229-643 



10 



15 20 
Log Frequency v (Hz) 

J0738-^1742 - PKS 0735-^17 



25 




4- 

4- 

























15 20 
Log Frequency v (Hz) 

J0746-^2549 - B2.2 0743-1-25 



10 



15 20 
Log Frequency V (Hz) 

J0818-f4222 - S4 0814-1-425 



25 



o 










o 

T 


















In 


4.^ 






OJ 

E 7 

"> 




* 

•"V 

I 


i* 

J 




f CvJ 

E V 
'S 




^ It 




Log 

16 -14 








Log 

16 -14 


/ ■ 







10 



15 

Log Frequency v (Hz) 



20 



25 



10 



15 20 
Log Frequency v (Hz) 



25 



Fig,28. The SED of lJy0537-286 (J0539-2839, top left), PKS 0548-322 (J0550-3216, top right), IRAS-L 06229-643 
(J0623-6436, middle left), PKS 0735+17 (J0738+1742, middle right), B2.2 0743+25 (J0746+2549, bottom left), and S4 0814+425 
(J0818+4222, bottom right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, 
Planck ERCSC and non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are 
shown in blue; literature or archival data are shown in light gray. 
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Fig. 29. The SED of OJ535 (J0824+5552, top left), 4C 71.07 (J0841+7053, top right), PKS 0851+202 (J0854+2006, middle left), 
B2 0912+29 (J0915+2933, middle right), S4 0917+44 (J0920+4441, bottom left), and PKS 0921-213 (J0923-2135, bottom right). 
Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck ERCSC and non- 
simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in blue; literature or 
archival data are shown in light gray. 
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Fig.30. The SED of 4C55.17 (J0957+5522, top left), IH 1013+498 (J1015+4926, top right), IRXSJ 105837.5+562816 
(J1058+5628, middle left), 4C01.28 (J1058+0133, middle right), PKS 1057-79 (J1058-8003, bottom left), and Mkn421 
(J 1 104+3812, bottom right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, 
Planck ERCSC and non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are 
shown in blue; literature or archival data are shown in light gray. 
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Fig.31. The SED of PKS 1 124-186 (J1127-1857, top left), PKS 1127-145 (J1130-1449, top right), B2 1128+31 (J1131+3114, 
middle left), S5 1133+704 (J1136+7009, middle right), PKS 1144-379 (J1147-3812, bottom left), and 4C 49.22 (J1153+4931, 
bottom right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck 
ERCSC and non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in 
blue; literature or archival data are shown in light gray. 
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Fig.32. The SED of 4C29.45 (J1159+2914, top left), ON 325 (J1217+3007, top right), PKS 1217+02 (J1220+0203, middle 
left), ON 231 (J1221+2813, middle right), PKS 1219+04 (J1222+0413, bottom left), and 3C273 (J1229+0203, bottom right). 
Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck ERCSC and non- 
simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in blue; literature or 
archival data are shown in light gray. 
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Fig. 33. The SED of PKS 1244-255 (J1246-2547, top left), PG 1246+586 (J1248+5820, top right). 3C279 (J1256-0547, middle 
left), Uy 1302-102 (J1305-1033, middle right), Uy 1308+326 (J1310+3220, bottom left), and GB6B 1347+0955 (J1350+0940, 
bottom right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck 
ERCSC and non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in 
blue; literature or archival data are shown in light gray. 
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Fig. 34. The SED of IWGAJ 1407.5-2700 (J1407-2701, top left) and of 3C 298.0 (J1419+0628, top right), CS0643 
(J1423+5055, middle left), of PG 1424+240 (J1427+2348, middle right), IRXSJ 145603.4+504825 (J1456+5048, bottom left), 
and PKS 1502+106 (J1504+1029, bottom right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data inte- 
grated over 2 months, Planck ERCSC and non-simultaneous ground based observations are shown in green; Fermi data integrated 
over 27 months are shown in blue; literature or archival data are shown in light gray. 
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Fig. 35. The SED of BZQJ 1507+0415 (J1507+0415, top left), 4C -05.64 (J1510-0543, top right), APLib (J1517-2422, middle 
left), PG 1553+113 (J1555+1111, middle right), WE 1601 + 16W3 (J1603+1553, bottom left), and OS -237.8 (J1625-2527, bottom 
right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck ERCSC and 
non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in blue; literature 
or archival data are shown in light gray. 
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Fig.36. The SED of 4C38.41 (J1635+3808, top left), NRA0512 (J1640+3946, top right), 3C345 (J1642+3948, middle left), 
MknSOl (J1653+3945, middle right), ARP 102B (J1719+4858, bottom left), and lES 1741 + 196 (J1743+1935, bottom right). 
Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck ERCSC and non- 
simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in blue; literature or 
archival data are shown in light gray. 
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Fig. 37. The SED of OT081 (J1751+0939, top left), S5 1803+784 (J1800+782, top right), PKSB 1830-210 (J1833-2103, middle 
left), PKS 1833-77 (J1840-7709, middle right), 2E 1908.2-201 (J1911-2006, bottom left), and PMNJ 1923-2104 (J1923-2104, 
bottom right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck 
ERCSC and non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in 
blue; literature or archival data are shown in light gray. 
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Fig. 38. The SED of OV-236 (J1924-2914, top left), lES 1959+650 (J1959+6508, top right), IJy 2005-489 (J2009-4849, middle 
left), PKS 2052-47 (J2056-4714, middle right), IJy 2126-158 (J2129-1538, bottom left), and S3 2141 + 17 (J2143+1743, bottom 
right). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck ERCSC and 
non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in blue; literature 
or archival data are shown in light gray. 
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Fig. 39. The SED of IJy 2144+092 (J2147+0929, top left), 4C 06.69 (J2148+0657, top right), PKS 2149-307 (J2151-3027, mid- 
dle left), BLLac (J2202+4216, middle right), 4C 31.63 (J2203+3145, botton left), and PKS 2204-54 (J2207-5346, botton right). 
Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck ERCSC and non- 
simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in blue; literature or 
archival data are shown in light gray. 
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J2209-4710 - NGC7213 J2229-0832 - PKS 2227-08 
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Fig. 40. The SED of NGC 7213 (J2209-4710, top left), PKS 2227-08 (J2229-0832, top right), PKS 2227-399 (J2230-3942, mid- 
dle left), 4C 11.69 (J2232+ 1143, middle right), 3C 454.3 (J2253+ 1608, bottom left), and PKS 2300-18 (J2303-1841, bottomright). 
Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, Planck ERCSC and non- 
simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are shown in blue; literature or 
archival data are shown in light gray. 
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J2327+0940 - PKS 2325+093 J2333-2343 - PKS 2331 -240 
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Fig.41. The SED of PKS 2325+093 (J2327+0940, top left), PKS 2331-240 (J2333-2343, top right), and lES 2344+514 
(J2347+5142, middle). Simultaneous data are shown in red; quasi-simultaneous data, i.e. Fermi data integrated over 2 months, 
Planck ERCSC and non-simultaneous ground based observations are shown in green; Fermi data integrated over 27 months are 
shown in blue; literature or archival data are shown in light gray. 
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Table 7. IRAM data. 









Flux density 


Flux density 








86.2 GHz 


142.3 GHz 


J2000.0 name 


Source name 


Obs. date 


(Jy) 


(Jy) 


J0217+0144 


PKS 0215+015 


2010-02-01 


1.91 ±0.18 


2.0 ± 0.4 


J0237+2848 


4C 28.07 


2010-02-01 


1.491 ±0.017 


1.27 ±0.08 


J0238+1636 


PKS 0235+164 


2010-02-01 


1.87 ± 0.13 


1.92 ±0.27 


J0238+1636 


PKS 0235+164 


2010-07-29 


1.46 ±0.22 


1.6 + 0.4 


J0238+1636 


PKS 0235+164 


2010-08-02 


1.45 ±0.16 


1.49 + 0.27 


J0238+1636 


PKS 0235+164 


2010-08-24 


1.39 ±0.13 


1.33 + 0.13 


J0433+0521 


3C 120 


2010-03-01 


1.864 ±0.016 


1.034 ± 0.001 


J0530+1331 


PKS 0528+134 


2009-08-27 


1.58 ±0.10 


1.22 ±0.10 


J0530+1331 


PKS 0528+134 


2010-03-22 


0.81 ±0.14 


0.52 ± 0.09 


J0530+1331 


PKS 0528+134 


2010-09-20 


0.88 ±0.12 


0.57 ± 0.06 


J08 18+4222 


S4 08 14+425 


2010-10-22 


0.84 ± 0.10 


0.94 ± 0.25 


J0854+2006 


PKS 085 1+202 


2010-03-22 


6.0 ± 0.6 


5.0 ± 0.5 


J0854+2006 


PKS 0851+202 


2010-04-26 


3.95 ± 0.17 


3.12 ±0.29 


Jl 104+3812 


Mkn421 


2010-04-26 


0.6 ± 0.3 


0.38 ± 0.05 


J1104+3812 


Mkn421 


2010-05-31 




0.50 ± 0.06 


J1256-0547 


3C279 


2010-01-31 


13.7 ± 0.6 


12,5 ± 1.8 


J1256-0547 


3C279 


2010-06-14 


16.6 ± 1.1 


12.3 ± 1.4 


J 1256-0547 


3C279 


2010-06-17 


13.4 ±0.8 


10.4 ± 1.1 


J 1256-0547 


3C279 


2010-07-01 


13.9 ± 1.1 


11.2± 1.6 


J 1256-0547 


3C279 


2010-07-05 


14.44 ± 0.22 


12.0 ± 0.5 


J 1504+ 1029 


PKS 1502+106 


2010-08-02 


0.7 ± 0.3 


0.49 ±0.17 


J1504+1029 


PKS 1502+106 


2010-08-24 


0.9 ± 0.4 


0.53 ±0.10 


J1635+3808 


4C 38.41 


2010-03-22 


3.31 ±0.23 


2.7 ± 0.4 


J1635+3808 


4C 38.41 


2010-08-02 


3.46 ± 0.27 


2.9 ± 0.4 


J1635+3808 


4C 38.41 


2010-08-24 


3.58 ±0.18 


3.05 ± 0.28 


J1642+3948 


3C345 


2010-03-22 


5.04 ± 0.20 


3.9 ± 0.4 


J1653+3945 


Mkn501 


2009-08-27 




0.51 ±0.09 


J1653+3945 


Mkn501 


2010-03-22 


0.70 + 0.14 


0.56 ±0.10 


J1653+3945 


Mkn501 


2010-04-28 


0.78 ±0.15 


0.579 ± 0.004 


J1653+3945 


Mkn501 


2010-05-31 


0.81 ± 0.07 


0.666 ± 0.026 


J1653+3945 


MknSOl 


2010-07-01 


0.94 ± 0.29 


0.93 ± 0.04 


J1653+3945 


Mkn501 


2010-08-02 


0.87 ±0.19 


0.72 ±0.18 


J1653+3945 


Mkn501 


2010-08-24 


1.00 ± 0.25 


0.80 ±0.16 


J1653+3945 


Mkn501 


2010-10-23 


0.789 ± 0.019 


0.746 ± 0.018 


J2202+4216 


BLLac 


2010-05-31 


5.0 ± 0.6 


4.2 ± 0.8 


J2202+4216 


BLLac 


2010-07-01 


6.5 ± 1.0 


7.9 ± 1.8 


J2202+4216 


BLLac 


2010-06-08 


5.58 ±0.19 


5.8 ± 0.5 


J2202+4216 


BLLac 


2010-06-17 


5.7 ± 0.5 


5.6 ± 0.7 


J2202+4216 


BLLac 


2010-06-21 


6.8 ± 1.0 


7.7 ± 1.2 


J2253+1608 


3C 454.3 


2009-11-21 


31.3 ±0.6 


32.4 ± 1.3 


J2253+1608 


3C 454.3 


2010-05-31 


31.4 ±0.9 


27.3 ± 2.6 


J2253+1608 


3C 454.3 


2010-06-08 


25.66 ± 0.24 


19.7 ± 0.4 


J2253+1608 


3C 454.3 


2010-06-17 


26.0 ± 1.4 


22.2 ± 2.5 


J2253+1608 


3C 454.3 


2010-06-21 


28.8 ± 1.1 


28.2 ± 2.2 


J2253+1608 


3C 454.3 


2010-07-01 


25.11 ±0.07 


26.3 ± 0.9 
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p. Giommi et al.: Simultaneous Planck, Swift, and Fermi observations of blazars 
Table 16. Fermi-LA^ data at one bin/decade simultaneous with the Planck observations. 



Flux density (10~'"phcnr- s'' McV~' J 



J2000.0 name 


Source name 


lOOMev- IGeV 


1 - lOGeV 


10-lOOGeV 


0017+8135 


S5 0014+813 


<0.33 


<3.0 X 10-^ 




0035+5950 


lES 0033+595 




<0.32 


<5.5 X 10-2 


0048+3157 


Mkn348 


<0.49 






0120-2701 


IJy 01 18-272 


<0.70 


(1.0+0.4) X 10-2 




0204-1701 


PKS 0202-17 


0.48±0.21 






0210-5101 


PKS 0208-512 


<0.64 


<1.3 X 10-3 




0214+5144 


GB6J 0214+5 145 


<4.8 






0217+7349 


IJy 0212+735 








0217+0144 


PKS 0215+015 


0.91+0.31 


<7.7 X 10-^ 




0221+3556 


IJy 0218+357 


<1.4 


<1.1 X 10-2 




0237+2848 


4C 28.07 


1.5±0.5 


<1.1 X 10-2 




0238+1636 


PKS 0235+164 


<1.9 


<7.6 X 10-3 




0336+3218 


NRAO 140 


<0.95 


<1.5 X 10-2 


<4.5 X lO-'' 


0423-0120 


PKS 0420-01 


1.5±0.5 


<1.1 X 10-2 




0433+0521 


3C 120 


<1.3 






0539-2839 


IJy 0537-286 


<0.59 






0738+1742 


PKS 0735+17 


<1.1 


(5.2+3.1) X 10-3 




0746+2549 


B2.2 0743+25 




<0.46 




0818+4222 


S4 08 14+425 


<1.2 


<1.4x 10-2 




0824+5552 


OJ535 


<0.39 






0841+7053 


4C 71.07 


<0.43 






0854+2006 


PKS 0851+202 


<0.81 


<1.1 X 10-2 




0915+2933 


B2 0912+29 


<0.46 


<7.8 X 10-3 




0923-2135 


PKS 0921-213 


<0.21 






0957+5522 


4C 55.17 


0.89+0.27 


(1.7+0.5) X 10-3 


(1.1+1.1) X 10-4 


1058+5628 


IRXSJ 105837.5+562816 


0.51±0.19 


(6.2±2.3) X 10-3 




1058+0133 


4C 01.28 


<0.71 


<7.5 X 10-3 




1058-8003 


PKS 1057-79 


<0.65 


<7.4 X 10-3 


<9.7 X 10-5 


1131+3114 


B2 1128+31 


<0.24 






1136+7009 


S5 1133+704 


<0.77 


<8.0 X 10-3 




1153+4931 


4C 49.22 


<0.48 






1159+2914 


4C 29.45 


17.8±0.3 


(2.0±0.5) X 10-2 




1220+0203 


PKS 1217+02 


<0.50 






1221+2813 


ON 231 


<0.65 


<5.7 X 10-3 




1222+0413 


PKS 1219+04 


1.2±0.4 


<7.3 X 10-3 




1246-2547 


PKS 1244-255 


<1.8 


<1.7 X 10-2 




1305-1033 


IJy 1302-102 


<0.55 






1310+3220 


IJy 1308+326 


<0.31 


(6.5±2.2) X 10-3 




1350+0940 


GB6B 1347+0955 


<0.43 






1419+0628 


3C 298.0 


<0.66 


<7.2 X 10-3 




1507+0415 


BZQJ 1507+0415 


<0.45 






1510-0543 


4C -05.64 


(45.5±7.6) X 10-2 


(21.3+2.9) X lO-'' 


<5.8 X 10-^ 


1517-2422 


APLib 


1.22±0.41 


<9.5 X 10-3 


<4.8 X 10-"* 


1555+1111 


PG 1553+113 


<1.1 


(14.6+3.3) X 10-3 


<7.8 X 10-5 


1625-2527 


OS -237.8 


<2.9 


<1.4x 10-2 




1640+3946 


NRAO 5 12 


<9.8 


<0.67 


(5.7+1.7) X 10-2 


1719+4858 


ARP 102B 








1743+1935 


lES 1741 + 196 


<0.31 




<2.5 X 10-2 


1840-7709 


PKS 1833-77 


<0.59 


<3.3 X 10-3 




1911-2006 


2E 1908.2-201 


<1.4 


<2.4 X 10-2 




1923-2104 


PMNJ 1923-2104 


1.97±0.67 


(8.78±5.4) X 10-3 




1924-2914 


OV-236 


<1.6 


(12.7±6.9) X 10-3 




2129-1538 


IJy 2126-158 




<6.1 X 10-3 


<4.2 X lO-'' 


2147+0929 


IJy 2144+092 


1.47±0.45 


<7.6 X 10-3 




2151-3027 


PKS 2149-307 


<1.2 






2203+3145 


4C 31.63 


<0.32 






2207-5346 


PKS 2204-54 


<0.61 






Column 1 and 2: source name; Column 3-9 


Fermi-LAI flux densities in three energy bands. 
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p. Giommi et al.: Simultaneous Planck, Swift, and Fermi observations of blazars 
Table 16. Fermi-LA^ data at one bin/decade simultaneous with the Planck observations. Continued. 



JZUUU.U name 


Source name 


Flux density (10 "'phcm -s 


' MeV ' ) 


100 Mev- 1 GeV 


1 - 10 GeV 


L\j — lUUOeV 




NnC 791 


<0.62 






2229-0832 


PKS 2227-08 


<0.81 


<9.7 X 10-3 




2230-3942 


PKS 2227-399 


<0.31 




<9.1 X 10-5 


2303-1842 


PKS 2300-18 


<0.45 


<6.9 X 10-3 




2327+0940 


PKS 2325+093 


<0.97 


<4.7 X 10-3 




2333-2343 


PKS 2331-240 


<0.83 






2347+5142 


lES 2344+514 




<9.9 X 10-3 





Column 1 and 2: source name; Column 3-9: Fermj-LAT flux densities in three energy bands. 
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p. Giommi et al.: Simultaneous Planck, Swift, and Fermi observations of blazars 
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p. Giommi et al.: Simultaneous Planck, Swift, and Fermi observations of blazars 
Table 18. Fermi-LA^ data at one bin/decade integrated over 2 months around the Planck observations. 



19000 nnmp 


^ n ff m m 
oUUH^C llalliC 


Flux density (10 ^"phcm ^ s 


Mev ) 


lOOMev- IGeV 


1 - lOGeV 


10 1 00 npv 

i u — i UU vJC V 






<0.41 


<3.8 X 10^3 


< / .U X lU 


r\r\A Qi-ii 


MKn j4o 


<0.19 


<4.6 X lO"'' 


Owl n— 5 
<0.o X lU 


VlZU—Z /Ul 


ijy Ul lo—ziz 


(43.3±9.8) X 10-2 


(4.7+1.2) X 10-3 


/"/i A_i_'j ^\ V/ 1 n— 5 
(O.U+J.j) X lU 


AO 1/1 1 C 1 /I /I 

Uz 14+3 144 


r^Yi t^f AO 1/1 1 C 1 /I c 

UdoJ (Jz14+M4j 


<0.22 


<1.9 X 10^3 




049^—0 1 90 


PK'*^ 0490—01 


1.34+0.18 


(6.5+1.5) X 10-3 




0738+1742 


PKS 0735+17 


0.51+0.13 


(3.6+1.2) X 10-3 


(7.3+4.2) X 10-5 


0854+2006 


PKS 0851+202 


0.50+0.14 


(3.2+1.19 X 10-3 


(4.1+2.9) X 10-5 


0915+2933 


B2 0912+29 


<0.18 


(96.5+5.4) X 10-5 




1305-1033 


IJy 1302-102 


<0.29 






1743+1935 


lES 1741+196 


<0.39 


<2.9 X 10-3 




2347+5142 


lES 2344+514 


<0.41 


<2.5 X 10-3 





Column 1 and 2: source name; Column 3-9: Fermi-hKV flux densities in three energy bands. 
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p. Giommi et al.: Simultaneous Planck, Swift, and Fermi observations of blazars 
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Table 20. Fermi-LA^ data at one bin/decade integrated over 27 months. 



Flux density (10-'"ph cm'^ s'' MeVQ 



J2000.0 name 


Source name 


100 Mev - 1 GeV 


1 - 10 GeV 


10- 100 GeV 


0017+8135 


S5 0014+813 


<8.7 X 10"^ 


<1.8 X 10 




0214+5144 


GB6J 0214+5 145 


<6.7 X 10"^ 


<4.2 X lO-'' 


<6.4 X 10-^ 


0217+7349 


IJy 0212+735 


(28.1+3.8) X 10-2 


(4.8+1.7) X 10-"^ 


<3.7 X 10"* 


0550-3216 


PKS 0548-322 


<5.8 X 10^2 


<4.7 X lO""* 


(4.6+2.8) X lO"*" 


1147-3812 


PKS 1144-379 


(23.1+3.0) X 10^2 


(15.7±2.2) X lO-'' 


<8.3 X 10"'' 


1220+0203 


PKS 1217+02 


(12.3±3.6) X 10-2 


<4.9 X lO""* 


<5.8 X 10-* 


1305-1033 


IJy 1302-102 


<6.0 X 10"2 






1419+0628 


3C 298.0 


<0.10 


<4.0 X lO-'' 




1423+5055 


BZQJ 1423+5055 


<4.2 X 10-2 






1456+5048 


IRXSJ 145603.4+504825 


<0.12 




<6.5 X 10 


1510-0543 


4C -05.64 


(45.5±7.6) X 10-2 


(21.3+2.9) X 10-"* 


<5.8 X 10"* 


1603+1554 


WE1601 + 16W3 




<2.5 X 10""* 


<4.9 X 10"* 


1719+4858 


ARP102B 


<4.8 X 10-2 


<1.4x lO-'' 




1743+1935 


lES 1741 + 196 


<7.6 X 10-2 


(9.5+1.9) X lO-'' 


(5.8±3.2) X 10-* 


2129-1538 


IJy 2126-158 


<0.11 


<3.5 X lO-'' 




2303-1842 


PKS 2300-18 




<2.7 X lO-'' 


<2.3 X 10-* 


2333-2343 


PKS 2331-240 


<9.5 X 10-2 


<4.3 X lO-"* 


<3.7 X 10-* 



Column 1 and 2: source name; Column 3-9: Fermi-'LPS flux densities in three energy bands. 
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